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Chapter 1 


Introduction 


1.1 Posing 


The purpose of this thesis is to contribute to the development of the mathematical 
theory of evolving large population. In view of its numerous applications — in life and 
social sciences in particular — this theory has become popular within the last decades. 
Due to the intrinsic character of the processes that take place therein, as well as 
due to the large size, the evolution of such populations is naturally considered as 
stochastic. This feature predetermines the ways and the means of the mathematical 
modeling of such objects. In the present thesis, we work in the Markov approach in 
which the population states are probability measures defined on appropriate phase 
spaces, the evolution of which is obtained by solving (in one or another way) the 
corresponding Kolmogorov-Fokker-Planck equations [6]. The first key aspect of the 
thesis is that the age of the population members (time of their presence in the 
population) is explicitly taken into account. The second aspect is that the studied 
populations are infinite, which opens the possibility to clearly distinguish between 
the local and the global aspects of the theory, see [25]. The latter peculiarity of the 
theory leads, however, to essential technical complications as compared to the case 
of finite populations. 

In view of the complexity of the evolution of large populations, its description is 
conducted in different spatio-temporal scales, cf. |2, 5]. At the macroscopic scale, the 
population is characterized by aggregated parameters, like density, mobility, etc. For 
such parameters, one derives — rather heuristically deduces — corresponding evolution 
equations intended to describe the population dynamics. The particular form of 
such equations corresponds to the model being considered. Below in this section, we 
outline some of them with the focus on the models containing the age parameter. 
The microscopic description is characterized by the use of the so-called individual 
based models in which each single population member is assigned its individual 
traits. In our case, such traits include also the member age. The most comprehensive 
description here amounts to constructing stochastic Markov processes. In certain 
sense, we work in the space between these two approaches. 


Most of the material of this thesis was published in my articles [18, 19, 20]. In 
view of this, certain parts of the present text are direct quotations of the correspond- 
ing parts of these articles, supplied sometimes with suitable comments. For the sake 
of reader’s convenience, some important formulas may be repeated. The majority of 
the proofs follow directly after the corresponding statements. The present thesis has 
the following structure. In the remaining part of this section we provide introduc- 
tory material and then outline the aims of the thesis and its main results. The next 
chapter is dedicated to providing necessary technicalities including the Markovian 
terminology as well as our way of describing member’s ages as marks. In Chapter 3 
based on the results of [18], we describe the evolution of a model introduced herein 
in terms of correlation functions. This might be considered as an intermediate ap- 
proach connecting the macro- and the microscopic theories. In Chapters 4 and 5 
based on [19, 20], we provide the microscopic description of the evolution of another 
model of this kind introduced by us. This model is somewhat simpler than that 
studied in Chapter 3. In this case, we directly construct the Markov evolution of 
states. In Chapter 4, we consider the case where the population habitat is X = IR^, 
which allows us to introduce a special class of the states — in fact, it is the same 
class as in Chapter 3, consisting of states possessing correlation functions — and then 
to show that the evolution leaves this class invariant. This can be considered as an 
additional information concerning the properties of the Markov evolution of this 
model. In Chapter 5, we study the case where X is just a locally compact Polish 
space. Here we prove the existence of a Markov process corresponding to this model. 
This is performed in the framework of the martingale approach [11, 13] combined 
with the theory of the Fokker-Planck equation [6] in the spirit of [26]. 


1.2  Age-structured population models 


The use of population models traces back to 1798 when one of the most known 
model of this kind was proposed by Thomas Robert Malthus. The author con- 
sidered a homogeneous population living isolated in an unchanging habitat with 
unlimited resources, in which the speeds of both procreation and mortality are pro- 
portional to the population size. This model predicts either an unlimited growth or 
an inevitable asymptotic extinction of the population. Due to these far from realistic 
assumptions and predictions it received rather restricted applications to real-world 
objects. Afterwards, a number of improvements and refinements were invented and 
implemented. Among such improvements we remark those which take into account 
the age structure of the populations being studied. Let us mention some motivations 
for this. In populations of living beings, individuals with different ages may have 
different survival capacities, mortality and reproduction rates. Many models of the 
infectious disease initiation and transmission include structuring the host population 
by the disease state dependent on the time since the disease initiation. Along with 
biological, biomedical and ecological applications, the age aspects are being taken 


into account also in many sociological and financial models, see [1, 4, 7, 21, 25, 29]. 
In particular, this kind of modeling is used to handling network security problems. 
W. Murray was the first who suggested to link modeling the ecology of the computer 
viruses to its counterpart dealing with the biological ones [29]. J. Kephart and S. 
White proposed to employ the SIS-kind models for studying the spread of computer 
viruses [21]. 

As mentioned above, the age structured populations can be studied at both 
microscopic or macroscopic scales. Below we outline some typical models of this 
kind representing each of these types. 


1.2.1 Macroscopic phenomenology 
The Lotka-McKendrick's model 


The Lotka-McKendrick's model is a direct analogue of the Malthus model. The 
population dwells in an invariant habitat and its members differ in age, see [15]. 
The fertility 8(a) and mortality u(a) rates are intrinsic parameters. In the simplest 
case, they do not depend on time. The evolution of the age-dependent population 
density p(a,t) is described by the following equation 


apla, t) + 5,p(a,t) = —n(a)pCa,t) 
p(0,t) = Jy B(s)p(s, t)ds (1.1) 
p(a,0) — po(a) 


where 2-p(a, t) is to capture aging — a uniform drift with unit speed along the age 
axes |0,a;] C R+. The basic and detailed analysis of this model can be found in 
[15]. This model laid the foundations for the forthcoming works in this field based 
in part on the theory of Volterra integral equations, for more details see [14]. 

Iannelli in [15] considered also some modifications of the Lotka-McKendrick's 
model, obtained by adding a time dependence of the vital rates G(a,t) and u(a,t). 
He studied the asymptotic behavior of these models, including strong and weak 
ergodicity, see Definition 3.1, page 37 and Definition 3.4 page 39 ibid. In Theorem 
5.4.4 below, we will also consider ergodicity problem for the Markov process in one 
of our models. 


The SIR model with aging 


In this version of the SIR model (Susceptible, Infectious, Recovered) of infection 
spread an isolated population dwells in an invariant habitat, structured by age. 
Similarly as in (1.1) p(a,t) is an age-density function at time t with age a € [0, a4), 
with possible a, = --oo. The population is divided into three groups: suscepti- 
ble, infected and removed, which are described by their respective age-dependent 
densities s(a, t), (a, t), r(a,t) at time t, i.e., 


p(a,t) = s(a,t) + i(a,t) + r(a, t). 


The evolution equation now reads 
Ss(a,t) + Zs(a,t) = —(u(a)  A(a,t))s(a, t), 


Zila, t)+ Žila, t) 


Ala, t)s(a, t) — (ula) + o(a)  9(a))i(a, t), 


rla, t) + 2-r(a, t) 


—p(a)r(a, t) + o(a)i(a, t). 


A detailed description of this model is provided in section 9.5 of [8], assuming that 
coefficients describe HIV infection and named it by analogy the SIA model from 
susceptible individuals, infective individuals, and AIDS cases. Then it was addi- 
tionally assumed that there is a minimum age of interaction between the members 
of the population. In Section 9.5.3 of [8], this model was studied by Co-semigroup 
methods. We will also use to these methods in the Theorem 5.3.8 below. 


The COVID modeling 


The very recent examples of the age-structured modeling are devoted to the world- 
wide spread of the COVID-19. Unfortunately, we have the opportunity to see this in 
the real life. During the COVID-19 pandemic, researchers were using and formulat- 
ing mathematical models as a technique in gaining insight into the mode of spread 
of the virus, e.g. [1, 33]. The authors of [1] investigated the situation where vac- 
cines and drugs could not be applied. Furthermore, they considered three different 
response strategies: total lock-down, partial lock-down, aiming at achieving herd 
immunity. They enriched the model SIR. with economic components capturing the 
outbreak of the pandemic. This led them to the SIR. model with wealth dynamics. 
The population was divided into two groups under and over the age of 65, assuming 
that members of the elderly group are more susceptible to infection. This model 
is very accurate and requires familiarization with many parameters, if the reader 
would like to delve into this topic, we recommend reading the entire article [1]. 


1.2.2 Microscopic individual-based theory: finite systems 
The Jagers-Klebaner model 


The authors of [17] proposed a modification of the known Bellman-Harris branching 
process model [3]. Therein, the particles are born and die at random. In the event 
of death of an individual, there arises a random number of independent offsprings. 
The birth parameters might depend on the population size and age structure. The 


generator G of the process is given by the formula: 


F'((F, AFA) + X balay HECO) + (0.49) - FC, A))} 


j=l 


at 2. hala; AEA[F(Y (az) f(0) + (f, A) — f(a5))] (1.2) 


GF((f, A) 


F((f, A))}, 


where F((f, 4)) is a test function on the space of measures u, f is a function defined 
on R, A = (a1ı,a2,...,az) is the finite collection of individuals ages, h(a) is the 
hazard rate, b4(a) is the intensity of birth, Y (a) the number of offspring distributed 
depending on the age structure of the whole collection .A and on the age of the 
particle which died, E4 the expectation number of individuals depending on the 
initial collection. The states of the model are finite positive Borel measures on 
R+ equipped with the weak topology. In [17], the authors used some martingale's 
techniques to establish the asymptotic properties of this model. 


The Méléard-Tran model 


In [28], Méléard and Tran introduced an age-structured individual-based model, 
similar to (1.2). They considered a finite population in continuous time t where the 
individuals reproduce and die with rates depending on a quantitative trait x € R and 
on their age a € R,. The trait x is a special feature that characterizes individuals, 
it can mean body size, rate of food intake, etc. The model is described by the 
generator L having the form 


LF(u) = (S f CF (u) + i 


XxR+ 


(d(x, a) 


+ pU(x,0)) (F(u — Sea) — Fit) 


+ b(æ,a)(1 — p) (Filu 9) — Fr) 


Weadp f (Fru + desno) — Fi) kle, hah] de, da) 


where X is a closed subset of R%, u is a finite measure on X, f is a continuous 
bounded real-valued function with bounded continuous derivatives with respect to 
a, d(x,a) is the natural death rate, F is a cylindrical function. A large and finite 
population is considered and the population size is represented by the parameter 
n. Furthermore the birth and death rates are taken of order n. In this model, the 
trait can change depending on the occurrence of mutations with probability p. The 
individual compete with each other, which is captured in the rate U(z, a). 


1.2.3 Microscopic individual-based theory: infinite systems 


The necessity of using individual-based models of infinite systems for modeling large 
real-world objects of this kind was realized in the first half of the XX-th century. It 
was then materialized in the concept of the ‘infinite-volume limit’ used in the statis- 
tical mechanics of physical particle system. In the present context, it was recognized 
and discussed in, e.g., |22, 25]. Briefly, the basic idea of this approach consist in 
getting the possibility to clearly distinguish between local and global aspects of the 
theory. As a finite system occupies a compact subset of the habitat (assuming the 
latter be locally compact), its description is always local rather than global. In deal- 
ing with infinite systems, one employs probability measures on an appropriate phase 
space to model system's states. Such models of states are quite abstract objects, 
not appropriate for a direct investigation. In view of this, it is convenient to restrict 
the theory to a class of states that possess correlation functions |24, 27], and then to 
describe the evolution of such states as the evolution of these correlation functions. 
The corresponding approach was realized in [5, 22, 25], see also the works quoted 
in these publications. Its additional advantage is that it is suitable for connecting 
micro-and macroscopic descriptions, as well as to find out intermediate versions of 
such theories. The only disadvantage is that the habitat should have some additional 
properties, e.g., it should be IR?. 


1.3 The aims and overview 


In this work, we introduce and study two individual-based models, described by their 
Kolmogorov operators, similarly as in (1.2), of an infinite population the members 
of which are characterized by their age - time since appearance in the population. In 
Chapters 3 and 4, as the habitat we take R? and use probability measures possessing 
correlation functions, see Definition 3.1.1. In this case, we describe these states with 
the help of such functions. In Chapter 3, based on [18], we introduce and study 
the model (3.1), where new members are born by the existing population members. 
This model turns to be quite hard to study. In view of this, we restrict ourselves to 
solving the evolution equation for the first two correlation functions in subsect. 3.3.1 
and 3.3.2, which are usually studied in various macroscopic theories. However, in 
our approach the evolution equations are directly deduced from an individual-based 
model. These results are formulated in Proposition 3.2.1. 

In Chapter 4, we introduce another model, see (4.1), in which particles arrive 
to and depart from X independently of each other. Their traits are (x,a), x € X 
and a > 0 being time since appearance. In this chapter, we take X = R? and 
use the same class of states as in Chapter 3. Due to the fact that the model is 
now simpler, we managed to prove in Theorem 4.2.2 that the corresponding Fokker- 
Planck equation has a global in time solution t — u+ which describes the evolution 
of states of this model. Additionally I found a stationary state for this evolution. 
This chapter is based on [19]. 


In Chapter 5, based on [20], I consider the same model as in (4.1) with the habitat 
X being just a locally compact Polish space. In this case, one cannot use methods 
based on correlation functions. Therefore, the very construction of the theory needs 
to be modified, including introducing special metrics on the corresponding configu- 
ration spaces. The latter proved to be quite technical. The main result is contained 
in Theorem 5.4.4, where the existence of a unique Markov process describing the 
evolution of this model was proved by showing that the corresponding martingale 
problem is well-posed. Additionally, assuming that the departure rate is separated 
away from zero I proved that this evolution is temporarily ergodic. 

To summarize: the main results of this thesis are 


1. Proposition 3.2.1 and the solutions obtained in subsects. 3.3.1 and 3.3.2. 
2. Theorem 4.2.2 and its proof. 


3. Theorem 5.4.4 and its proof. 


Chapter 2 


Preliminaries 


2.1 The Markov evolution 


2.1.1 Generalities 


Throughout this work we use the following notations: N = {1,2,...}, No = NU {0}; 
R7. d € N, is a standard Euclidean space. By a Polish space we mean a separable 
topological space the topology of which is consistent with a complete metric. Mostly, 
such a space is locally compact, which means that each of its points has a compact 
neighborhood. For a suitable A C E, by 14 we denote the indicator of A, i.e., the 


functions such that 
1 ifzrzecA4, 
la (x) := i (2.1) 


0 if rgs. 


For a Polish space E, by B(E) we denote the corresponding Borel o-field; C,(£) 
(resp. B,(£)) stands for the set of all bounded and continuous (resp. bounded 
and measurable) functions f : E — R. By C} (E) and Bj (E) we denote the set 
of positive elements of C,(E) and B,(E), respectively. Finally, C4(E) C (E) 
consists of all continuous compactly supported functions. 


Definition 2.1.1. A family of functions, F, is said to separate the points of a Polish 
space E if for each distinct x,y € E, one finds f € F such that f(x) 4 f(y). 


By oF we denote the smallest sub-field of B(E) such that each f € F is oF- 
measurable; by P(E) we mean the set of all probability Borel measures on (E, B(E)). 
For z € E, the Dirac measure 6, with center at this x is defined by its values 


&(A)-1aA(2,  A€B(E), (2.2) 


see (2.1). For a given measure u and a suitable function f, we write u(f) = f fdp. 
Then 6,(f) = f(x). For a sequence {Un}nen C P(E), by writing u, > u € P(E) 
we mean its weak convergence, i.e., Un(f) > u(f) for all f € (E). 


Definition 2.1.2. A family of functions F is said to be separating if ui(f) = pa( f) 
holding for all f € F implies uj = uo for any uj, ua € P(E). 


If F separates the points of E and its linear span is an algebra with respect to 
pointwise operations, then it is separating, see Ethier and Kurtz (1986), Theorem 
4.5 on page 113, [13]. 


Proposition 2.1.3 (Ethier and Kurtz). Let V and F be a complete and separable 
metric spaces and a family of functions F : V — IR, respectively. Assume that: 


(a) each F € F is bounded and continuous; 
(b) Fi, Fy € F, their pointwise product is in F; 
(c) for each distinct v4, v9 € V , there exists F € F such that F(vi) # F(v3); 
(d) F contains F = 1. 
Then F is separating. 


Definition 2.1.4. A family of functions F is said to be convergence determining if 
Un( f) — u( f) holding for all f € F implies un = u as n — +00. 


2.1.2 The Kolmogorov-Fokker-Planck formalism 


Let E be a Polish space. Broadly speaking, a dynamical system on F is defined by a 
family of maps {T;}:50, Ty : E — E, such that T; o Ts = T;,, holding for all t, s > 0, 
and Tọ = I. Here o and J denote composition and the identity map, respectively. 
The points x € E are then treated as states of the system; i.e., x, :— Txo is its 
state at time t. In the stochastic (Markov) version, states are probability measures 
u € P(E), and thus u = Tuo. Then the ‘point’ states r € E appear here as 
dx € P(E), see (2.2). That is, the evolution of the system is said to be deterministic 
if the family (71,50 preserves the set of point states, which means that T;ó, = ôs, 
for all t > 0 and x € E. 

In the Kolmogorov-Fokker-Planck formalism, see |6, 22, 23], the construction of 
the evolution u — pz is performed by solving the Fokker-Planck equation 


i F) = i (F) + f m(LF)s, m € P(E), (2.3) 


where L is the so called Kolmogorov operator specific for the model being consid- 
ered. It acts on the test functions, also called observables, usually taken from a 
suitable subset of C (Æ). In fact, the choice of the domain of L predetermines the 
quality of the description of the stochastic evolution of a given model. Sometimes, 
one approaches to constructing the evolution in question indirectly by solving the 
Kolmogorov equation 


d 
qi I LF,, Filio EE Fo, (2.4) 


10 


which might be considered and ‘dual’ to (2.3). Then solutions of (2.3) are obtained 
from the solutions of (2.4) by setting u(F) = u(F;), where u and F are taken as 
the initial conditions of (2.3) and (2.4), respectively. Finally, a Markov process, 
X, with state space E is obtained by constructing a probability measure, P, on a 
space of paths with values in E such that the corresponding probability of the event 
X(t) € A, A € B(E), is given by (A), where ju is a solution of (2.3). This most 
comprehensive description will be obtained in Chapter 5 below by solving the so 
called martingale problem for L. 


2.2 Measures on configuration spaces 


2.2.1 The configuration spaces 


In the description of the stochastic evolution of infinite populations they are viewed 
as random ‘clouds’ of point ‘particles’ placed in a suitable habitat. In our case, 
the latter is a locally compact Polish space X. In Chapters 3 and 4, we take 
X = RA, whereas in Chapter 5 it is just a general locally compact Polish space. Each 
mentioned cloud is supposed to be locally finite, which means that its intersection 
with a compact A C X is supposed finite. The mathematical model of such a cloud 
is a configuration — a counting Radon measure on X. The set of all such measures I 
is then defined by the property that y(A) € No for each compact A C X. This can 
also be given the following geometric interpretation. Let 6 be a complete metric of 
X consistent with its topology. Then for each r > 0, one finds r’ € (0,r) such that 
the ball B? (x) = (y € X : (x,y) < r'} has compact closure. Therefore, for a given 
y ET and x € X, there exists r > 0 such that 4(B2(x)) < oo. Then we define 


p(y) = fw € X : inf(Bi)) = ny (a) € N}. 


Each y is fully characterized by the pair (p(y), n), in which p(y) is the ground 
configuration for y. It is convenient to extend n, to the whole X by setting n(x) = 0 
whenever x ¢ p(y). The mentioned geometric interpretation of y is the collection 
of the points x € p(y) such that at this x there is located n,(x) indistinguishable 
‘particles’ - population members. We refer the reader to [9, 10, 16, 27, 34] for more 
detail on this issue. 

Our aim now is to define a suitable topology on T. It is the weak-hash (vague) 
topology. By definition, it is the weakest topology that makes continuous all the 
maps 


ray | fend) - D nef Die), feCalX), (25) 


z€p(y) rey 


where the latter equality is the definition of 5... , See [27] for more detail. It is 
know, see [34], that [ equipped with this topology is a Polish space. Next, we define 


D-(.4€r: n(x) = 1 for all z € p(y)}. (2.6) 


11 


That is, y is in Dif y = p(y) as sets. Such configurations are called simple. The set 
T is a Gs-subset of T, and hence is a Polish space in the subspace topology, see [34]. 

Along with infinite configurations, I contains also finite ones. Let Io stand for 
the set of all such configurations. It can be viewed as a subset of the set of all finite 
Borel measures on X, closed in the weak topology. Note that the subspace topology 
induced on To be the weak-hash topology introduced in (2.5) coincides with the 
latter. At the same time, each y € I' can be interpreted as the element of X” /Xm, 
where m = |y] is its cardinality and Xm is the corresponding symmetric group. We 
repeatedly exploit this interpretation below. 

In the sequel, we will need the following extension of the last equality in (2.5). 
For y € I and z € p(y), by y \ x we understand the configuration y’ such that 
n (y) = n4(y) for all y Z x and: (a) p(y’) = p(y) and ny(x) =n, — 1 if n,(x) > 1; 
(b) p(y’) = p(y) Nx is n(x) = 1. Similarly we define y = y U y by setting: (a) 
py) = p(y) U ty}, ny (y) = 1 if y € p(y); (b) p(y’) = p(y) and ny(y) = n4(y) + 1 
if y € p(y). Thereafter, we have 


f(z,y) = f(z,y)y(dz)y(dy) — | f(z,x)y(dz), (2.7) 


that can also be generalized to all m € N 


yo oe y f (a, mg) (2.8) 


L1EY x2EY\21 LmEY\{21,---;Lm—1} 


| 
—— 

| 
rm 
— 
Q 


felur- , Yna Vd) +++ (dng) 


GEKm KE 
= ` (-1)' ` ae 25. felu, arini Yna), 
where Km is the collection of all graphs with vertices {1,2,...,m}, lg and ng are 
the number of edges and the connected components of G, respectively, whereas 
falyı, ---, Yng) is obtained from f(x1,..., £m) by setting the arguments 2),,...21,. 


of the latter equal y; where l4, . . ./,; are the vertices of the j-th connected component 
of G. Note that the summations in the left-hand sides of (2.7) and (2.8) have the 
direct meaning if y € Č, see (2.6), where y can be considered as just a locally finite 
subset of X. 

Now for y € l and m € N, by means of (2.8) we define the following counting 
measure on X” equipped with the usual product topology. For a compact subset 
A C X”, we set 


QUAS CS No caue ` lA(n,..., Em): (2.9) 


LIEY T2EY\T1 Sm YN sXm-1) 


12 


That is QU? (A) is the number of ordered strings (z1,...,r,,) contained in A. It 
is know, see [27], that the map y 5 QU? (A) is measurable for each compact A. 
However, it may be unbounded. 


2.2.2 Measures on Il 


Let B(T) be the Borel o-field of subsets of I relative to the weak-hash topology 
defined in (2.5). Recall that P(T) denotes the set of all probability measures on I'. 


Remark 2.2.1. As a Gs-subset, Ÿ is in B(T). Hence, each u € P(T) with the property 


i (T) = 1 can be redefined as a probability measure on the Polish space I. 

As mentioned above the maps y e QU" (A) need not be bounded hence p- 
integrable for a given u € P(T). Following [27] we will say that u has correlations 
of order m if 


(A) = f tym) < o, (2.10) 


holding for all compact A C X™. Then the Radon measure mu on X” is called the 


m-th order correlation measure corresponding to u. It is clear that the existence of a 
Y m/ <m. We say that u has correlation 
functions of all orders if it has all corresponding uU). 

Let us now consider the case of X = R¢. For m € N, let A? be the Lebesgue 


measure on X", i.e., 


given vU"? implies the existence of all wf" 


M™ (day,...,d%m) = dz: dtm. 


Assume that a given u € P(T) has mu for this m, which is A"-absolutely contin- 
uous and the corresponding Radon-Nikodym derivative ke is \")-almost every- 
where bounded. That is, Kh) is a symmetric element of L°(X™) := L~(X™, A0), 
Then we call ke the m-th order correlation function of u. By the very definition, 
it follows that 
k(x, erstma) = 0; 

The first two correlation functions have a special meaning characterizing the cor- 
responding state u. For a given compact A C X, let us consider OAJ, which 
counts the points of y contained in A, see (2.9) and (2.10). Then wP (A) is just the 
-expected value of the number of points contained in A. At the same time, 


w(D(A) = jJ TE KO (a)dz, 


which means that kP is just the particle density in state u. In a similar way, by 
(2.7) one shows that 


cule, y) := KP (x,y) — KL (a) ky? (y) (2.11) 
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is the truncated spatial correlation function in state u. These two functions kP and 
C, are usually employed in describing population dynamics at macro-level by kinetic 
equations. Similarly as in Lemma 2.10 of [26], one proves the following statement, 
in which I’ stands for the set of simple configurations defined in (2.6). 


Proposition 2.2.2. Let X = R? and u € P(T) have ke. Then p(T) = 1, ie., the 
configurations in state u are almost surely simple. 


The importance of this statement can be seen in the light of Remark 2.2.1. Let 
us outline the arguments yielding its proof. For e € (0, 1), one defines 


1 
|x — y| 


Hy=)) > h(ny, = Ay) = 


TEY yey\x 


It is clear that H(y) = +00 for each y € T \ Č. At the same time, by (2.9), (2.10) 
we have 


uan = | QP hud) = [Pehle y)dedy < oe 


The latter holds true as ke € L^? (X?) and h is absolutely integrable. It implies 
uT V D) = 0. By similar arguments one can also extend Proposition 2.2.2 to the 
case where X is just a Polish space. That is, wT) = 1 if u possesses the second 
order correlation measure verifying wh (D) = 0 where D = ((z,y) € X*: £ = y). 
Finally, let us mention a probability measure on I which plays a special role in 
this work. Let o be a positive Radon measure on X. The latter is just a general 
locally compact Polish space. The Poison measure 7, with intensity measure @ is 
defined by its correlation measures 
Qu es pm (2.12) 


which in particular means that for A = A; x Am, all A; C X compact, one has 
wl) (A) = o(Ai) -+ e(As). 


That is, in state 7, the points choose their positions in X independently of each 
other. If X = R? and o(dr) = p(x)dx, p € L®(X), then all ken exist. By (2.11) 
and (2.12), it follows that 


G2 y) = 0; 


i.e., the spatial correlations are absent in this state. 
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2.2.8 Configuration spaces with marks 


As mentioned above, our principal aim is to take into account also particles ages 
— the time of their presence in the population. To this end, we pass to the traits 
$ = (x,a), x € X anda € Ry := [0,+00). It is then naturally to consider 
configurations 4 as collections of such compound traits in which œ appears as a 
mark. More on the general theory of marked configuration spaces can be found in 
[10]. According to this, a marked configuration ^ is a counting Radon measure on 
X := X x R, such that, for each compact A C X, the set A = A x R} verifies 
4(A) € No. That is, every compact A contains only a finite number of the elements 
of 4 of different ages. 

In Chapter 5, we deal with the general case where X is just a locally compact 
Polish space, without assuming a priori properties of the states. This means that 
the corresponding marked configurations may be multiple, and thus each x € X is 
characterized by a finite (possibly empty) configuration of ages a = {a1,...,@}. In 
Chapters 3 and 4, we deal with populations dwelling in X — IR^ and with the states 
possessing all correlation functions. According to Proposition 2.2.2, such states are 
supported on simple configurations, which means that only those will be taken into 
account, see also Remark 2.2.1. In this case each age nonempty configuration a is 
a singleton. In view of this, we will denote there the age variable by a. In view of 
the mentioned differences in the properties of the marked configurations spaces, we 
make a more detailed presentation of them an the beginnings of the corresponding 
Chapters. 
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Chapter 3 


The birth-and-death dynamics 


In this chapter based on [18], we introduce a model and study the stochastic dynam- 
ics of an infinite population of point particles dwelling in X = IR?, which amounts to 
the following. The particle having trait $ = (x,a) at time t = 0: (a) can die at time 
t with probability 1 — e^" 9, where the death rate m is a suitable function; (b) can 
give birth to a new particle located in a compact A with probability 1 — e-": 9t. 
where v is a positive measure kernel on (B(X), X). As mentioned above, we restrict 
our consideration to states possessing all correlation functions and hence satisfying 
the conditions of Proposition 2.2.2. In view of this and for the further simplicity, in 
this chapter we denote by [ the space of all simple configurations defined in (2.6). 
Then the space of marked configurations D consists of the pairs 4 = (y,a), where 
a: y — IR, is the age map that assigns a, to each x € y. Its topology will be made 
precise below. 

The Kolmogorov operator corresponding to the model studied herein has the 
following form 


(FIG) = SFG) + Y mi)tFG Và) - FG (3.1) 


Here the first term on the right-hand side corresponds to aging, whereas the remain- 
ing ones describe the death and the birth of new particles as mentioned above. The 
model corresponding to (3.1) can also be used to describe the infection spread in an 
infinite population. In this interpretation, the population of healthy individuals is 
present as a vacuum, whereas ^ corresponds to the population of infected individu- 
als. The recovery depends only on the time since infection and is described by the 
function m(£). The spread of infection from the infected population is captured by 
the kernel v. 

For this model, assuming that the evolution preserves the set of states possessing 
all correlation functions from the Kolmogorov equation (2.4) we derive the evolution 
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equations for the correlation functions and solve them for k® and k®), which yields 
the result usually obtained from phenomenological kinetic-like equations. However, 
in our case we obtained it directly from the individual-based model described in 
(3.1). 


3.1 The technicalities 


As mentioned above, in this chapter the configuration space D defined as the set of 
pairs (y, a) where y is a simple configuration on X = R? and a: y —^ Rx. The value 
of a at a given x is denoted by a,. The space X = X x R, is equipped with the 
usual product topology and the Lebesgue measure d$ = dada. 

The configuration space P is equipped with the following topology. Let F be the 
set of all bounded and continuous maps f : X x R+ — R such that each f vanishes 
for x € A^ := XX A for an f-specific compact A C X. Then Î is equipped with the 
weakest topology that makes continuous all the maps 


P3460 M f(t, a2), JEF. 


rey 
With this topology Î is a Polish state. 


Definition 3.1.1. By P, we denote the subset of P(L) consisting of those u that have 
all correlation measures ae”, each of which is absolutely continuous with respect to 
the corresponding Lebesgue measure AC (dí, ---dêm) = dij d£, such that their 
Radon-Nikodym derivatives kem have the properties: for all m € N the following 


holds 
(i) ke is a symmetric element of L'e(X") = Lee(X, X6): 


(ii) kC? (4, 01) ..., (m, à)) — 0 as max; a; — +00; 


(i) for almost all zi,...,z4 € X and a1, ...am € R4, the map 


(0, +00) > t KO? (x1, a1)... , (Em, am + t)), 


is continuously differentiable; 


(iv) the following holds 


i KO) (Gian) vas (Op) darse didis E L®(X™); 


m 
(v) there exists a -dependent £ > 0 such that, for each 0 € L'(X) := L'(X, dê), 
|@|| = 1, the series 
oo gm "" | l l 
1+ D m! A K (24, ..., 2,.)0(23) ++ 0($5)d£i- - - d, 
m=1 


is absolutely convergent. 
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Now we introduce the following function 


Fa =] [aea yer. (3.2) 


yey 


Definition 3.1.2. By Oo we will denote the set of continuous functions 0 : X x 
R+ — (—1,0] having the following properties: (a) each 0(x,a) vanishes for x € 
A* for some 6-specific compact A C X; (b) for each x, the map a |> 0(x,a) is 
continuously differentiable. 


For 0 € Qo, is measurable and F^ satisfies 
0< F*(4) <1. 


Hence, it is u-integrable for all u € PÒ). By employing Proposition 2.1.3, it can be 
proved that the set of functions {F° : 0 € Qo} is separating. According to item (v) 
of Definition 3.1.1, for u € P, this function can be u-integrable also for 0 € L'(X) 
for sufficiently small ||@||. Indeed, by (2.8) and (2.10) it follows that 


u(F?) 


(ef [Xm XD 560962 aem (3.3) 


«(ze)» 


NEY LEN 


oo 1 "" l 
= OSEE ME: M4... $0 (21) «+ O( By) dB - dêm 
m-—1 


= Join oon 


rcr 


Here 7 € y means that 7) is a finite and n subset of 4, Îo is the subset of 
Î consisting of finite configurations, k,() = k (2,,...,4,,) for f) = {41,.-., 2m}, 
and A is the Lebesgue-Poisson measure on Îo defined D the formula 


[ GJA) = G(g) + +> as GO (2),, ..., 2, )d£i +++ Bm, 


where G and G/?? are defined analogously as k and k™ above. In the sequel, we 
repeatedly use the following statement known as Minlos' lemma [30]. 


Lemma 3.1.3. Let n € N,n 7 2, then for all measurable functions G : Pp > RA: 
To — R it is true that: 
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We also use the following evident formula 
2,4(8) 5, B= DD AGBA), (3.5) 
LEY ney\x ney Len 


holding for suitable functions h and G. 


3.2 The model 


To be able to guarantee that the evolution preserves the states possessing corre- 
lation functions with the properties as in Definition 3.1.1 we impose the following 
conditions on the kernel v that appears in (3.1). Recall that by a measure kernel we 
mean a map B(X) x X > (4,0) 5 vol X), which is a finite positive Borel measure v; 
for each § € X and a measurable function à ++ v4(A) for each A € B(X). Namely, 
we assume that, for each y, the measure v; is such that 


vy(de) = b(g|a)de, (3.6) 


where the map (ĝ, x) — b(g|v) € R+ is measurable and satisfies 


Va c R, b(y, alx) € B(x — y), f eae =: Oe: (3.7) 


Vz E€ X | b(g|x)dj € b < oo. (3.8) 
x 


The condition in (3.7) is to control the procreation of the new members uniformly as 
to the age of the parental particles. The second condition controls the procreation 
from all age groups to given compact subset of A C X, i.e., 


[Cf iioi) ae <m, 


where |A| is the Euclidean volume of A. 
Concerning the function (x,a) — m(x,a) we assume that x — m(x,a) is mea- 
surable for each fixed a and a++ m(x, a) continuous for each x. For 7 € l'o, we then 


set 
M(H) = 3 mê). 


Our first aim is to pass from the Fokker-Planck equation (2.3) to the evolution 
equation for the correlation functions. Here one may follow the way elaborated in 
[26]. For k : T'o — R, we then define 


(LK) (À D (0) - M (4)&(8). 
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In order to use this operator on a regular basis, one would have to introduce an 
appropriate Sobolev-like Banach space of such functions keeping in mind that they 
satisfy the conditions mentioned in Definition 3.1.1. In this chapter, however, we 
restrict ourselves to deriving only equations for k“) and k), which are usually 
‘deduced’ by employing heuristic arguments. Keeping this in mind, we prove the 
following statement. 


Proposition 3.2.1. Assume that u € P, is such that its correlation function k, 
satisfies the renewal condition, cf. [15], 


ks eu (0) = | ble wA\ EU GAG + kylA\ 8) Y Gln). — (3.9) 


yen 


Assume further that both integrals 


I0 := u(LEF®’) (3.10) 


absolutely converge for all 0 € O. Then I,(0) = J, (0). 


Proof. Define 


(L^k,) (f) 


(L^k,)(0) +X Slaa) (3.11) 


xrEN 


x 


kyla) + f bless Neo idi + bul \ 2) D (Glo | 


yen\x 


where d(a,) is Dirac’s ó-function, which means that the right-hand side of (3.11) is 
a distribution. Note however, that L^k = L^k if k satisfies (3.9). The Kolmogorov 
operator (3.1) with v as in (3.6) acts on F^ (3.2) as follows 


(LFf(4) = 5 [] 0+ 6@) +35 maa- aea [[ à (9) 


rey yer rey yer 


+ | Sai (ay)(1 + 6(9) — 1) [ [0 - e(5)43. 


LEY ZEY 
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Now we split LF? into three parts. The first one is 


furom = [OS™ TE a eonun 


LEY i 


- [X M Teer 


Ge ce yen 


- [3X2 II «6» 


Ay 


ncy xen yen 
= [sey eS DT eonun) 
Eo vey ^ yenw 


" [ (f mau KM 


= -f n2 (N U 2)0(a4) ul U ĉ)0 (2)d$ [ [6(9)^ 
To 


yen 


- -j ( 35995, M oaas) T] (9). 


0 xen yen 


This is obtained using (3.5) with A(x) = ND Bn) = Tye, 9(9) and Lemma 3.1.3. 
The second part is 


1 Ly F*(4)u(d¥) = f S "(6 &) J| a 9 e(9)u(a5) 
r rey E 
` -f Vom) S> D[otnucas) 
rey ncy\x YEN 
= -|E EmA Tua 


= x S m(#) (HH) T] OOA) 


zc yen 


21 


Finally, the third part 


framan = f f Eea [Ta + dina) 
r TEY ZEY 
3 (pe b(2|9)6(a,)0(G)(1 + 9(@)) [T (1 + (2) 9449) 
zey\x 
= ines b(z|y)d( \(1 + e($ D SS Tf e@dguca4) 
ncy\x 2€n 
= rre (ay)(1 +0(ĉ2 ) lH 0(2)düu(d^) 
ncy xen zen\xUy 
= [fm ñ) X blêg) (ay) I] 6(2)dyX(dn) 
Do JX ren z€nNUy 
NEL "TONS (ay) D[ AA ñ) 
LEN zenuy 
= So )b(#|9)5(ay) | [96a2X(48) 
To yen zen 
+ f SERAN )b(&|y) 0 (ay MICE )dyX(dn). 
lo yen zen 


Putting everything together and taking into account (3.9) we get 


J (Er tat tat nut = [ ie (ñ) (IIoc )) aa 


ren 
= | L^k,(fj e (IIoc )) ac, 
To ren 
which yields the proof, see (3.10). 
3.3 The result 
As in [26], by Proposition 3.2.1 one may obtain solutions of the Fokker-Planck 
equation (2.3) by solving the Cauchy problem 
o AG. act 
zh) = (LR), kdo = E, (3.12) 


Ot 
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where k,,, is the correlation function of the initial condition of (2.3). As mentioned 
above, a regular theory of this equation should be based on the use of appropriate 
Banach spaces. After developing such a theory, one should also prove that k; is 
indeed the correlation function for a certain state, which is nael a hard technical 
problem, see [25]. Here we obtain from (3.12) the equations for k; ^ and ki which 
can be viewed as the usual PDE of transport type. 

For T > 0, set 


v. (1)) = kee) (3.13) 
where 
th = {(@,d¢ — T) : (£, az) EH}, T < min{az, tj. 


To get a classical solution of (3.12) we have to assume that k+ is continuously differ- 
entiable with respect to each a, and therefore it is continuously differentiable with 
respect to t. By taking the derivative in (3.13) by (3.12) we get 


d . A " 
which can be solved 
elf) = vui yexp f MGydi). 


By (3.13) this yields, where a; — min Qx: 
zen 


ko(1],) exp (- {iM (Ña d), t <a, 


kta, (flag) exp (— fo? M(fjg)d9) , t > as, 


ki) = (3.14) 


which is a functional equation. Below we study this equation in the two particular 
cases mentioned above. 


3.3.1 The evolution of densities 


By setting 7) = ($) in (3.14), we get the following formula describing the evolution 
of the first correlation function 


kK (x, a — t) ) exp (- - - 6)d0) , t<a 


Gas (3.15) 


KC (x,0) exp (— fom x ' t a, 


subject to the following renewal condition obtained from (3.9) 


kK (2,0) = f b(gle)k (9)dj. 
X 
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By (3.15) we have 


kf? (a, 0) 


f (| bly, a, |x)kf? (y, ay)da) dy 
R4 0 


t 
| f Msc, (0, e mms qaa (3.16) 
R4 Jo 


+00 
+ f f b(y, ay |) ks? (y, ay — tye Jo mav- da, quj. 
R4 Jt 


By (3.7), the second term in (3.16) is bounded uniformly with respect to t and x by 
- ess sup | [ow 
ycR2 
see also item (iv) of Definition 3.1.1. Set 


KY (a, 0) := ulz). 


Now (3.16) can be written in the following form: 


u = Au +v, (3.17) 
where 
(Au)(x) = [fe (y, ayz) Ut—ay (y je m(y,a, —0) 4 da, dy, 
+00 j 
u(r) = | / b(y, a, |a) kd (y, ay ie- Jo myay— da, dy. 
R4 Jt 


Both u and v are positive elements of the Banach space Ua of functions R, x IR? — R, 
which are C! (R+) for almost all fixed x and L*(I4?) for all a, equipped with the 
norm 


llulla = ED e?" lus ros qa); a > 0. (3.18) 


By (3.7), the operator norm of A satisfies 
|All < Besssup | e 09 Jo myb) da =: Bqla). 
yeR? JO 


If the mortality rate m(y, a) is separated away from zero, i.e., satisfies m(y, a) > Mx, 
then 
qla) € 1/(a + m,). (3.19) 


For Bq(a) < 1, we can solve (3.17) in the form 
uc(f-A)' 
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where v is determined by the initial condition of (3.12). By item (i) of Definition 
3.1.1 and (3.8), it follows that 

ella < [ello S ISP le co (3.20) 


Now (3.16) in the space Ua takes the form: 


ki (2,0) = w(x) = (Bv), (x) = D> (^v), (2). (3.21) 


n=0 


Finally, the evolution of the particle density is described by the formula 


ko (a, Gate Jo m(z,a—8)d0. t€ 

n ee . 

ki (zx, a) ` (Av) a (x)e~ Jo ena Ea (3. ) 
n=0 


The solution (3.22) may increase in time. By (3.18), (3.20) and (3.21), for t > a it 
satisfies 


b 
(za) € — e eat « FHI cet. 


"lBeta) ~ 1- 8qla) 
For m, > B, by (3.19) we have Gq(0) < 1, which means that kK (x, a) remains 
bounded in this case. 
3.3.2 The evolution of the second correlation function 


In this case, we take /j = {(x, az), (y, ay)} with az < ay. By (3.14) we then get 


x (JEW Jo (™(y,ay—8)-+m(x,a2—8))d0 tSas 
ke ay dag (3.23) 
ke (y, Ay — às, 2, 0)e^ Jo” (m(y,ay- 8) m(z,as—0)d0. t> ds. 


To find KO (y, à, — ax, x, 0), we use (3.9) which yields 


-Foo 
kO(y, a, v,0) = b(z, a; ax )kC? (y, a; z,a,)da,dz (3.24) 
' R4 Jo 


+ kK (y, a)b(y, ala). 


25 


For t < a, by (3.22) for kC (y, a) and (3.23) for k” (y, a; z, ae) we bring (3.24) to 
the following 


bly, al) KO (y, a — tje fo ma-od6 (3.25) 


+ / f tn yk 2 (y, a= az; 2,0) 
m4 Jo 


x exp ( — [ [m(z, az) + m(y, ay — 6)de) da;dz 


-Foo 
+ J / b(z, a; |a) KS (y, a —1;z,a; — t) 
R4 Jt 


x exp ( — A [m(y, a — 0) + m(z, az — 8)] dd) da;dz. 


ke (y, a, T, 0) 


On the other hand, for t > a, we rewrite (3.24) as follows: 


kP (y,a,£,0) = bly, alx) (Bv), , (y)e- lo me (3.26) 


ff (z, az|£) ki”, (y,a— a4; 2,0) 
Ra 


x exp (- E [n(z, a; — 8)) + m(y, a — ay da;dz 


t 
+ f f Melk usos isa, — a) 
R4 Ja 


jo eu (- | E Hea 2) dico 


+oo 
+ / f b(z, a-lx)k® (y, 0; z, az — a) 
R4 Jt 


x exp (-f [n(y,a — 0) + m(z, a, — oja) da,dz. 

0 
Now we would like to solve (3.25) and (3.26) as a single equations in the space of 
continuously differentiable functions R} 3 t œ> w, € C! (R,) @L®©((R2)?) @L1(R,). 
Let: 


wily, a, x) TT k,O (y, a, T, 0), 


then we have: 


wily, a, x) = (Aow),(y, a, x) + fily, a, x). (3.27) 
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For t < a define: 
$ 
C E = I I A E ae 
R4 JO 


x exp (- [ont a — 0)) + m(z, a; — oao) da,dz, 


+00 
E 1 n edo Gato (3.28) 
R4 Jt 


x 


t 
exp (-f m(y,a — 9) + m(z, a, — oio) dady 
0 


+ 


ky (y, a)b(x, y, a), 


and for t > a: 


(Aow).(y, a, x) = | f b(z, a;|z)wi s, (2, y, a I az) 
R4 JO 


exp (- [ (m(y, a — 0)) + m(z, a, — 0)d6) da;dz 


x 
+ f. [re az|r)Wr_alY, Z, a; — a) 
x exp ( — [mw a —0)) - m(z, az — 6))a0) da,dz 


0 
+ I f b(z, az|£)WialY, 2, az aa a) 
R2 Jt 


x exp ( — [on a —0))+ m(z,a, — 6))a0) da,dz, 


FAY a) = k, O (y, a)b(y, a, £). (3.29) 


Let W, be a space with the norm: 


lwll = sup e *' esssup_ | |[w:(y, a, x)|da. 
£20 —— @y)e(R4y? 70 


It is clear that: 


f - (Ant), a, )da = I oaa f AG saat 
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Therefore 


oo t a 
f (Aqw):(y, a, x)da = f il f b(z, 0, X), a, (z,y,a m az) 
0 0 JR4 J0 


— (-[ m(y,a — 8)) ^ m(z,a, = 2) da,dzda 
0 


t t 
+ 1 Um 
0 JRi Ja 


x exp (- / EUER yan) da,dzda 


t oo 
+ f 1 f b(z,az|£)Wi-aly, z, a; — a) 
0 JRd Jt 


x exp (- / Oia yan) da,dzda 


oo t 
+ f f f b(z,az|£)Wia,(2,Y,a — az) 
t R4 JO 
x exp (- [ m(y,a — 0)) — m(z, a; — yan) da,dzda. 
0 


To estimate the previous integral we use (3.8) and the fact that e ^ < 1, where 
a>0. Leta = A for some fixed q < 1 therefore: 


lwla < allele: 


We can write the solution for (3.27) as 
k(x, 0, y, a) = = w( T,y,0 =X ( (AQ” P T.y, a). 
n=0 


The solution for k? (x, as, y, ay) takes the form: 


k (x, a, — t; 2 — t) 


x exp ( E Áo (m (£, az — 0) + m(y, ay — 6)a0), for t < dy; 
(2) = eS 
k`“ (£, az, Y, ay) = b» (A f), a, (y, ay — az, £) + kD, (y, a)b(y, a,|x)) 
n=0 


x exp ( — fo" (m(x, a, — 0) + m(y, ay — 0)d0), Tor az < t. 
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For a, < t < a, with the appropriate formula for f, (3.28): 


Mafias) = ff Meals (esa a) 


jc (- i: Moa) nee Qr) dude 


For ay < a, < t, f, like in (3.29): 


(Afitnae) = ff Messi Geta — a) 


x 


T (- ^ Ge ieee nw) dude 


4 i f enis. — a) 


T (- f a E eae 22) dé 


" f f E EAT 
R4 Jt 


emm (- / emot adem esae: Qr) diss: 


0 


x 


x 
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Chapter 4 


The independent appearance 
model: evolution in a spacial class 
of states 


In this and the next chapters we introduce and study the following modification of 
the model defined in (3.1). The corresponding Kolmogorov operator is taken in the 
form 


(LF)5) = Yar) + m(é) [FN 4) - F())] (4.1) 


where the first two terms have the same meaning as the corresponding terms of L 
given in (3.1), whereas the third one described an independent appearance of new 
particles in the habitat governed by the measure x. In this chapter, we take X = IR" 
and 


x(dy) = b(y)dy, (4.2) 


with an appropriate positive density b and the Lebesgue measure dy. In the next 
chapter, we consider the case of a general locally compact Polish space X and hence 
a general arriving measure x. Our present choice of X and x (4.2) will allow us to 
solve the Fokker-Planck equation (2.3) and thus to construct the evolution of states 
Lo — Ht Which preserves a special class of them. We begin describing this class by 
introducing ‘tempered configurations’. 


4.1 Tempered configurations and measures 


In this section, we select a subset of D, C P such that the measures in question have 


^ 


the property L(D.) = 1, and thus forget of the remaining configurations. 
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4.1.1 Tempered configurations 


Let Y : X — R+ be: 


(a) continuous, bounded and strictly positive; 
(b) integrable, i.e., satisfy fy v(r)dx < oo. 


One can take y(x) = e~! as an example of such a function. Then we set 


wy=S ve), Fel, (4.3) 


which can take infinite values for some 4. Thereby, the set of tempered configurations 
is defined as : . 

D,—Í[y€r:w(4) « oo), 
where Î is the set of all (possibly multiple) configurations. Similarly as in [26], 
subsect. 2.3, we equip this set with the following metric 


v. (5,7) = sup Y sewe) -Y swa) 


REY REA! 


| (44) 


where the supremum is taken over the subset of the set of bounded Lipschitz- 
continuous functions C7" (X) consisting of those g : X — R for which 


= / 
sus SC «up lg (x, @) sty, o?) 2% 
E azgex It — y| - loc — a'l 


It is possible to prove, [26], Proposition 2.7, that the metric space (X, v,) is complete 
and separable. Let BÈ.) be the corresponding Borel o-field of such subsets of am 
By Kuratowski's theorem (Theorem 3.9 in [31]), one then proves that T, € B(Î) 
and BÈ.) coincides with the Borel o—field related to the topology on Ô, induced 
by the vague topology of D. This allows one to redefine each he PÈ) with the 
property u(Î,) = 1 as a measure on (Î,, B(P,)), see [26], Corollary 2.8 for further 
details. 

Let P, C P(Î) be as in Definition 3.1.1. By Proposition 2.2.2, we know that 
each u € P, is supported on the set of simple configurations, which allows us to 
assume — in this chapter also — that Î and Î, consist of single configurations only. 
It turns out that the set of measures has another useful property. 


Proposition 4.1.1. For each p € Py, it follows that u(P.) = 1. 


Proof. Similarly as in Proposition 2.2.2, the proof consist in showing that (WV) < oo, 
which readily follows by (2.9) and (2.10) 


u(t) = f. kP (é)b(a)dé < / w(a)dx (sup i ! nr) « oo, (4.5) 


X rcX 


see item (iv) of Definition 3.1.1. 
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. As mentioned above, each u € P, can be redefined as a probability measure on 
(L., B(V.)) with single configurations, which we assume from now on. Recall that 
the functions F are defined in (3.2). 


Definition 4.1.2. For a given measurable q : X — (0,1) and u € P(L), the measure 
u3 defined by the relation u*(F?) = u(F®) is called an independent q—thinning of 
u. Here 0,($) = 0(2)q(2) and 0 € Oo, see Definition 3.1.2. 


For u € P,, its q-thinning amounts to multiplying k,,(7) by Ib q(£). An im- 
portant subclass of P, constitute Poisson measures 7, with intensity measures such 
that 

o(d$) = p(2)d4, (4.6) 


with p satisfying the conditions of Definition 3.1.1 in the part related to kP. For 
such measures, one has, cf. (2.12), 


ks) = | [26 (4.7) 
Then by (3.3) it follows that 


To(F°) = exp ( f p(@)0(@)d2). (4.8) 


X 


Note that each u € P(P) can have a correlation function understood as a distribu- 
tion. To see this, let us first define 


; R vaslê —Goyy), if lf] = Êl = n; 
d ay ll [3 


0 otherwise. 
In the first line, £ — (21,...,24], f) (du... Gn}, En is the symmetric group, 


and o(2 — ĝ) is the usual Dirac ó—function on R^ x R}. The correlation function ky 
of the ó—measure à; € P(Î) is then 


k = X aÊ). (4.9) 


By (4.9) we then have 


&(F*) = | kô T] o(@) (a8) = [[a + ot. 
2 $cé 


BEY 


By means of k; one can define the correlation function for any u by the formula 


kl) = f BÒNA): (4.10) 
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4.1.2 Convolution of measures 


For u, u2 € P(L), their convolution is defined as 
(ave nay(P)= f f PGs oa eon 53) (4.11) 
in 


As mentioned above, the class F°, see (3.2), we 0 € O is separating hence measure- 
defining. Thus, it is enough to define ju; « u2 for such F°. By (3.2) it readily follows 
that 


(11 * ua) (F^) = p (F^)us (F^). (4.12) 
Then for u1, u2 € P, one readily gets that 
quu) — | (E OR D) OA. (418) 
To *écp(a) ven 
Proposition 4.1.3. It follows that p x u2 € Pa, whenever m, uo € Px. 
Proof. We begin by slightly formalizing Definition 3.1.1. According to its item (iv), 
it follows that each 


PAUCOS k®) ((2, a4), ..., (n, Gn))day -++ das, n EN, 
u dos d 
+ 


is essentially bounded. Let P ) be this bound. Then item (v) of Definition 3.1.1 is 
equivalent to the following 


nw « nle” 


n S : n EN, 


holding for some € > 0. At the same time, by (4.13) it follows that 


Rises apo ( DP ky (ñ \ Êk u2 (&), 


&Cp(1) 


which then yields 


xm < ME nw et le "xn < ZO 
m 
m=0 


H1*H2 pi puz 2 


where € = min[e,,;6,,). This yields the validity of item (v) of the mentioned 
definition. 


4.2 The result 


In this section, we formulate and prove a statement describing the evolution of our 
model introduced in (4.1), (4.2). It is determined by the model parameters m and b 
which are supposed nonnegative, measurable and bounded. Additionally, we assume 
that the map Ry 2 a — m(z, a) is continuous for each fixed zx. 
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4.2.1 The statement 


Along with the class of functions Og introduced in Definition 3.1.2, we will use the 
following class. 


Definition 4.2.1. By © we mean the collection of all 0 : X — R that have the 
following form 


O(a, a) = O(x)e- TPH) E e "V G)é(o) _ 1, (4.14) 


where V : X — (—1,0] is a continuous functions with compact support, i» is as in 
)- 


(4.3), T 2 0 and o(a 


Eres 


Then we define F = (F^:0 € ©} with F? having the form 


F'(4) = exp (3 7 los(1  9()) - 73 5 9(e(0)). 


Our aim is to solve the corresponding Fokker-Planck equation (2.3) for F? with 
0 € O. Note that 0 < F?(4) < 1 for each 4 € P, and u(F°) < 1 for all u € P,. It is 
possible to show, cf. [12], Theorem 18, that each F° is v,-continuous (see (4.4)). The 
pointwise product of F and F” is F% with 0" corresponding to 0" (x) = (a) * 9 (x) 
defined in (5.23) below and 7" = T +7’. Assume that 4, 442, both are in P,. Then 
one finds $ which belongs to exactly one of these configurations, say 4,. If there 
is no ĝ € 4» with p(y) = p(ĉ), one takes 7 = 0 and V such that 9(p(2)) Z 0 and 
Ü(p($)) = 0 for all € 42. Otherwise, one takes 7 > 0 and 9(p(2)) = 9(p($)) # 0 
and V(p(2)) = 0 for all 2 € 4; U 42 such that 9(p(2)) Z 0(p(£)). In both cases, the 
corresponding F° separates 44 and 4», see property (c) of Proposition 2.1.3. Clearly, 
F? = 1 for r = 0 and ? = 0. Then by Proposition 2.1.3 F = {F° : 0 € ©} is 
separating. 
Let us prove now that LF? is u—integrable for each u € P,. By (4.1) we have 


(LF*)\(4) = E (366.0) = mts 0)0(e,a)) PG \ 2) (4.15) 


+ F'G) | WAE, 0)dz =: H (f) + H(^). 


Since b is bounded, H» is also bounded. Since 9 is continuous and compactly sup- 
ported, it is -bounded. Hence, by (4.14) one concludes that, for all ¢ € X, the 
following holds 


MOIS wle), | (2), (4.16) 
where dg depends only on the choice of 0 and 7. By (4.16) we then have 


|71(7)| € DoW), (4.17) 
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holding with an appropriate Dy. By (4.5) this yields the property in question. Now 
for 0 € O and m as in (4.1), we set 


a+t 


0,(,a) = 0(x, a + t) exp ( — f m ajda), t>0, (4.18) 


Q 


and then define a map P, 3 u — ut € P,,t > 0 by the following relation 
(eo) = u( F”), 0 € 0o. (4.19) 


Since the family (F^ : 0 € Ol is separating, each ji is uniquely determined by 
(4.18), (4.19). Note that the correlation function of u^ can be expressed through 
that of u as follows 


ku (0) = Ji) ku) exp (-xf m m(z, a) da). (4.20) 


where ff = {(x,a — t) : x € p(f)), 


zu Ji(z,o0), Ji(z,o) := 1— 1po(0), 


By Definition 3.1.1, the map u p preserves P, and is a combination of a thinning 
and an age shift. Now we are at a position to formulate our result. 


Theorem 4.2.2. For each ug € P}, the Fokker-Planck equation (2.3) has a solution 
of the following form 
Ht = Ho * Tor, (4.21) 


where ut, is obtained from uo according to (4.19) and T, is the Poisson measure, 
see (4.6), (4.7) and (4.8), with the intensity measure 


o. (dx, da) = b(x, a)1)0,4)(a)drda := b(x) exp (- i m(x, 9«) 1j, (o) dxda, 
0 


(4.22) 
where 1jo2)(a) is defined in (2.1). If m($) > m, for some m, > 0, the evolution 
given in (4.21) has a stationary state T, with o(dx,da) = b(x,a)dxda, see (4.22). 


4.2.2 Comments 


Let us make some comments to this statement. According to (4.12), (4.8) and (4.21) 
it follows that 


(F°) = exp (aa b(r, o) 8r, a)dedor) wo F*). (4.23) 


39 


Hence, the solution satisfies the initial condition j4|i2o = Ho, see (4.18). If uio( 2) = 1, 
i.e., the initial state is an empty habitat, by (4.23) it follows that ui = m, with o; 
given in (4.22). Let us show that this u, satisfies (2.3). For 7,, by the following easy 
to prove formula, cf. (3.4), 


[ 57 G(En\ Ade) = » l G(E,n) (ae) (AA), 


9 cn 


we have that 


To(Hı) 


n (IIo) 2 [s - m(a,o)|&() l[ ^» — (24 


zen LEN yen 


= [ (IIo) fE moo] oco ET oov 


zen yen 


sí " p(z, 0)8(s, 0)d) (F^) 


— (fa [E + mle, a) od) P) 
And also 
To(H2) = ( if b(z)6(r, 0)d Jm. (F^). (4.25) 


In the sense of distributions, we have that 


o ð 
8a Loo (0) = — Ato): 


Then for p;,(x, a) = b(z, a)l (a), see (4.22), one obtains 


E +m(a, o) | pula, a) = - S. o). (4.26) 


By (4.24), (4.25) and the latter equality it follows that 


Tal LEF’) 


exp ( P p 80d) © [ AOON (4.27) 


= 2 e( f aaas), 


by which one readily concludes that ju; = To, satisfies (2.3). 
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4.2.3 The proof 


The proof of Theorem 4.2.2 is divided into the following steps 


(i) proving that for each 0 € O, the map t 4 ju(F?) has a continuous derivative 
at each t > 0 (by Lebesgue's dominated convergence theorem); 


(ii) proving that this derivative satisfies (4.28), see below; 

(ii) showing that 7, is a stationary state; 

(iv) proving the weak convergence of probability measures on je 
(v) 


The proof of the first part will be done by showing that: 


showing that the family (ut)t>o is tight (by Prohorov’s theorem). 


(a) for each 0 € O, the map t > ju(F?) has a continuous derivative at each t > 0; 


(b) this derivative satisfies, cf. (4.27), 
a G7) = we LP"). a2) 
By (4.11), (4.12), (4.21), (4.22) and (4.23) we have 
ILE?) = pol F® Jra (F°) =: olEF” Qilt). (4.29) 


In view of (4.27), the continuous differentiability in question will thus follow by the 
same property of t ++ uo(F®). By (4.18) we have 


ZF^4) = E (Ge8)r*GV5) E me, a&G)F^6 V4 
ve€p() ze€p(3) 
= M a(£)F^(5 8) =: (5). (4.30) 
ve€p(4) 


Similarly as in (4.17) we then conclude that ope «)| < D,V(¥), with a certain 


D, > 0. By Lebesgue's dominated convergence theorem this yields 


Sono FP) = n (F^) = (8) (4.31) 


f E (o. TT. e)a, 


zep(f) yep( Ni) 


| 
=> 
a 
£e 
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as well as the continuity of the map t +> po( ZF ^). Here kp is the correlation 
function of jug understood in the sense of (4.10). Now let us turn to proving (4.28). 
By (4.29) and (4.31) we have 


LHS(4.28) = uo(S;)Qe(t) + mE") f b(x)0(x, t) exp ( — i ma, <)ds) de. (4.32) 


X 


At the same time, by (4.15) it follows that 
Hi (1 U 42) = Hin) F^ (52) + Hy (âo) F’ 9), 


which by (4.11) and (4.21) yields 


RHS(3.17) = uh (H1 )To (F°) + p (F^) (H3) + pC F’) i b(z))(r)de, — (4.33) 


Note that 
Ta (F5) = exp ( L pos) = Qut) (4.34) 


see (4.8), (4.22) and (4.29). By (4.15) we have that 


= M hs) (2)F*(5N 2) 


TEY 


hg(z, a) := a0. a) — m(x, a)0(z, o). 


By (4.20) one then gets 


Mn) = f tuens) TT o) (4.35) 


ren yen 


[Cf scio moa) [Tono enn 


yen 


= [ff moe 


x exp ( — ] mea) ho(x, a)deda) 


AÀ) J [2u ev) exp ( - f 


yen Ay—t 


ay 


x 


my, <)ds J he(y, a)deda) 


I 
A 
T 
=> 
wn" 
> 
Ze 
e 
cass. 

M 
Lm 

a 
3 
wn" 
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Here ĝt and J; are as in (4.20) and 0, is defined in (4.18), whereas 


Ed e [ / yg lf U (7.0 — hair, a) exp ( — ] mins casae 


[ou 


-" f [feto tute o (- f m(:r,<)ds) drda 


By (4.18) and (4.30) we have 


as m(xz, 9«) = g(x, a). 


ho(aya+thexp(— | 


Q 


We use this in the latter expression and then in (4.35); thus, we arrive at the 
following 


usan) - f (IT IT AJA) 75). — (4.36) 
To LEN yen\x 


see (4.31). Now similarly as in (4.24) we obtain 


mH) = ( Í: holê)oldâ)) Ta (CF) (4.37) 


= (f, [ oo (7 f mesa) 
ð 


x [5 —m(z, a)| O(a, adda) mo, GF?) 


= (Jo | exp ( — f me os). n = (x, 0) | dx) m, (F^). 


Finally, we use (4.36) and (4.37) in (4.33), take into account (4.34) and (4.32), and 
conclude that (4.28) holds true. 

To prove that 7, is a stationary solution of (2.3) we again use (4.15) and (4.24). 
For 


b(x, a) = b(x) exp ( — A m(x, b)dð), 
we have, cf. (4.26), 


E 4 mz, e) b(x,a) — 0, 


which by (4.24) yields 7,(LF°) = 0, and hence the property in question. 
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Chapter 5 


The independent appearance 
model: a Markov process 


In this chapter we study the model defined in (4.1) with X being just a locally 
compact Polish space. In this case, we do not have the property as in Proposition 
2.2.2. That is why, we will deal with spaces of multiple configurations and then 
understand the corresponding sums as in (2.8). Here we also will not use tempered 
configurations and thus special classes of measures. The presentation of the material 
is close to our paper [20] and the notations used here are mostly independent of the 
previous part of the thesis. 

Since the configurations here are multiple, we take this into account and instead 
of (4.1) we write the following formula for the corresponding Kolmogorov operator 


(LF)O = 3 aO $5 mz, 0) [F4\(@,0))- F()]) (5.1) 


(x,a)€4 (z,a)e4 
" f [F(4 U (2,0)) — F(4)] x (da). 


The model parameters are subject to the following assumptions: 


(i) The departure rate X x R4 2 (x,a) > m(z,o) € Ry is continuous and 
bounded, i.e., such that m(x, a) < m, for some m, > 0 and all (xz, o). More- 
over, there exits « : [0, 1] — R+ such that &(c) — 0 as e > 0 and the following 
holds 


Veex|m(x,a) — m(z,o^)| € «(|a — a'l), la — o/| € (0, 1). (5.2) 


(ii) The arriving measure x is just a positive Radon measure. 


The result of the present chapter can be outlined as follows. We introduce a Banach 
space C of bounded continuous functions F : T — R, in which we define L as a 
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closed and densely defined linear operator that satisfies the conditions of the Hille- 
Yosida theorem, and hence is the generator of a Co—semigroup (S(t)hzo. Then 
the solution of (2.4) is obtained in the form F, = S(t)Fs. For a class of functions 
Fo, Fi corresponding to Fo € Fo is obtained in an explicit way. This allows us 
to explicitly construct the corresponding Markov transition function p} and obtain 
finite-dimensional laws of a Markov process X with values in D, which describes 
the stochastic evolution of our model. Possible objects of this kind are specified as 
stochastic processes that solve the martingale problem for L. Then we show that this 
problem is well-posed, i.e., uniqueness holds. The main ingredient of the proof here is 
showing that the corresponding Fokker-Planck equation for L has a unique solution, 
which we do by employing the resolvent of L. Assuming that m(x, a) > mo > 0, we 
also show that the process X has a unique stationary state, explicitly constructed 
in the paper, such that the laws of V(t) weakly converge to this state as t + +00. 


5.1 The space of marks and the metric 


The space X is equipped with the product topology assuming that the topology of 
R, be defined by the metric which we introduce now. 
For a > 0, we set 


w(a) = min { 3 
and then 


r(o,0) = wla), (5.3) 


r(a,ao) = min {la — oo[;w(a) + w(ao)|), 


where || is the usual absolute value of 6 € R. 


Proposition 5.1.1. The above introduced r is a metric such that (IR,,r) be a com- 
pact metric space. 


Proof. 'To prove the first part we just have to check the validity of the triangle 
inequality 

r(01,02) € r(o1, 03) + r(as, 03). (5.4) 
This technical exercise is made in Appendix. To prove the compactness, we have to 
show that: 


(a) r is complete; 


(b) the space (IR, r) is totally bounded. 


Let {a@n}nen C R+ be an r-Cauchy sequence. 
Here one may have the following possibilities: 


(i) there exists a < oo such that a, € à, such that n € N; 
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(ii) the considered sequence contains a subsequence that diverges in the usual 
sense. 


In case (i), {Qn }nen contains a subsequence, say {Anz }nen, such that [ang —a.| — 0 
as k — +00 for some a, € à. At the same time, for € < 2. r(O4, 04) < € implies 
[An — Am| < £, see (5.3), which means that r(an,a,.) 4 0 as n > +00. 

In case (ii), the divergent subsequence converges in r to zero, which implies that 
the whole sequence converges to zero in r. Hence, the latter metric is complete. To 
prove (b), we set B-(a) = {œ € R+ : r(a, a’) < e}. Fixe € (0,1) and take the least 
k € N such that k +1 > i. Then R, — Ls B-(je), which yields the property in 
question. 


Let us now compare r with the absolute-value metrics of R+. By C(IR,, 7.) we 
will mean the set of all bounded r-continuous functions, whereas C, (IR.,, 7j) is going 
to stand for the set of all bounded | - |-continuous functions. 


Proposition 5.1.2. 7, is coarser than T, and hence the embedding (R4, Ta) — 
(R+, 7.) is continuous, whereas both latter topological spaces are Borel isomorphic. 
Moreover, 


CR, 7.) = {u € lR, Ta) : Em u(a) = u(0)). (5.5) 
Proof. The validity of the first statement and (5.5) readily follows by the fact that 
each |- |-convergent sequence is also r-convergent, and each r-convergent sequence 
either converges in |- | to the same limit a ¥ 0, or has two |- |-accumulation points: 
0 and oo. Since the mentioned embedding is continuous and injective, it is also 
Borel-measurable. By Kuratowski's theorem, see Parthasarathy (1967), page 21, 
[31], its inverse is also measurable and thus is the isomorphism in question. This, 
in particular, means that the corresponding Borel o-fields coincide. 


Definition 5.1.3. For a suitable u € C(IR,, T.), we introduce the map a — u'(a), 
meaning the usual derivative for a >0 and the right-hand side one if a = 0. 
A given u is said to be continuously differentiable on R} if u' € C(IR,, T,). 


Let us consider the following functions 


2 


aeER,, necN. (5.6) 


eae | 


It is clear that: 


(a) each u,, is continuously differentiable, see Definition 5.1.3; 
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VA 


2 
b) un is d ing f > 4/—.M , unla) < d 
(b) un is decreasing for a E oreover, unla) < To an 
2a — nat c 
: = —— Más : == 5.7 
Un(@) (tna)? Jun (a)] < Vw (9.7) 


the latter holding for some c > 0 and all a > 0. 
Now let {0} }ken =: X C [0, +00) be countable and such that: 
(1) 01 = 0; 
(ii) Veen Ok < Ok41; On — 0 < œ as n — +00. 
Next, for k,n € N, we set 
Wn pla) = e THU), (5.8) 
Then wkn is continuously differentiable and the following holds 
lvi e(0)] = rta), (5.9) 
where C is the same as in (5.7). 
Next, let a be a finite collection of points o; € R}. That is, a = {a1 }i<i<m, 01 € ota 
for all |. For a € a, by nala) € N we will denote the multiplicity of o in a, i.e., the 
number of elements of a coinciding with this a. We extend it to all a > 0 by setting 


nala) = 0 whenever a is not in a. Two such a and a’ are equal if they consist of 
exactly the same elements, with the same multiplicities. 


Proposition 5.1.4. Let a anda’ be as just described. Then they are equal if 


do wrala) = Y 7 na ON (5.10) 


aca aca! 
holding for all k,n € N. 


Proof. For a as above and ¢ € C, consider 
faal) = »3 e $n) n EN, 
aca 


Each such f is an exponential type entire function. By (5.8) and (5.10) we have 

that (fna — fa, )|x = 0, holding for all n € N. Since X has a limiting point, this 

implies f, (6) = fa," (C) for all Ç € R and n € N. Obviously, Qn CE NN 
—+00 


where na(0) > 0 is the multiplicity of a = 0 in a. Then the just mentioned equality 
yields na(0) = na (0) and also 


pe ee ee (5.11) 


a€aM(01 o€a/N(0) 
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Let a, and a’, be the least positive elements of a and a’, respectively. Take n > =, 
a_:= min(a,; o/ ). Then, for such n and all a > a_, one has un(a) > un(a). Now 
we multiply both sides of (5.11) by e6“‘? and pass to the limit ( — —oo. This yields 
that o, = oz and nalas) = Nna (a). Thereafter, we subtract the coinciding terms 
from both sides of (5.11) and proceed to comparing the remaining least elements of 
a and a’. This eventually yields the equality to be proved. 


Let A be the set of all a = (ou]ixiza; m € No, 0 € o1 € a2 € ... Om. 
Define 


pa) = Y, 7 ra) (513) 


knéN 1+ pin(a, a’) 
Pea (a.a) =| wrala) — Y s wanla)| 
aca aca’ 


By Proposition 5.1.4 it follows that p is a metric on A. Each a € A can be 
considered as a finite counting measure defined on the compact space (R,,1r), for 
which a(A) = J aca Lala) = |a U A], holding for all Borel subsets A. The weak 
topology of A is defined as the coarsest topology that makes continuous all the maps 
a — J aca Wla), w € C(R,,7,.). In the weak topology, A is a closed subset of the 
space of all finite positive measures on C(R,, 7,). 


Proposition 5.1.5. (.A, p) is a complete metric space. The corresponding metric 
topology coincides with the weak topology that turns A into a locally compact Polish 
space. 


Proof. As each wz, is in C(R,, 7.), the weak convergence of a sequence (a4, sew C 
A to a certain a € A yields p(a, am) — 0,m —> +00. Assume now that {am}men is 
a p-Cauchy sequence. By taking o = 0 we then get from the latter that, for some 
m, € N, the cardinalities of all am, m > m,, coincide. By Prohorov’s theorem this 
yields that {am }men contains a subsequence that weakly converges to some a. Hence, 
the whole sequence converges in p to this a. Then the metric is complete and the 
corresponding metric topology is exactly the weak topology of .A. The separability 
and local compactness follow by the fact that C(IR,, 7.) is compact. 


Let T° be the set of all locally finite simple configurations on X. That is, each 
y €T is a subset of X such that each compact A C X contains a finite number of 
the elements of y. Let now ^ be the pair (y,n),y € T and n : — N. The value 
of n at a given x € y can be considered as the multiplicity of x € ^. That is, 7 
is a configuration with multiple locations, for which y is the ground configuration. 
Sometimes, we will write n4(x) to explicitly indicate that we mean the multiplicity 
of x in the mentioned 4. By Č we denote the set of all such multiple configurations. 
For 4 = (y,n), we write y = p(¥). The weak-hash, vague, topology of Č is defined 
as the coarsest topology that makes continuous all the maps 


y2 Mi nla)g(z),  sg€C«X) 


z€p(3) 
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It is well-known, see [34], Lemma 1.2, that with this topology Č is a Polish space, 
whereas Î is a Gs subset of I’, by which it is also Polish. Following Lenard [27] we 
will also consider ^; as configurations of point particles, in which distinct particles 
may have the same location. Such particles can be enumerated, which allows one to 


write 
3 n(x)g(x) = M g(a), (5.13) 
z€p(3) LEY 


where in the second sum we mean a certain enumeration of this sort. In the same 
sense, we will write 


3 ` mor ` g(xi,...,2,), MEN, 


LEY r€Nroi LEY\{L1,£2,.--;Lm—1} 


where in expressions like 7 V x we treat x the singleton {x}, whereas ¥ \ x is the 
measure such that na ,(r) = n3(r) — 1 and ns\2(y) = ns(y) for y A x. 

It is known, see, e.g., Zessin (1983), page 397, that there exists a collection 
{vs tsen =: V CC CZ(X) of suitable functions such that the metric 


44)-Y LED 4G) =|Lae@)- Law), — (629 


is complete and consistent with the weak-hash topology of I’. In the sequel, we will 
always mean this topology of n Obviously, we can and will assume that V contains 
also the following functions. 

Let 6 be a complete metric of X and X’ a countable dense subset of X. Each 
x’ € X' has a countable base of compact neighborhoods, which we denote by D(z’). 
Each A € D(z’) contains balls 


BAG) = {x E€ X : (x, x) <q} 


with compact closures, where q is a rational number satisfying q < q' for a ó-specific 
q € Q. For z/ € X',6 € D(x’), q € qd ands € (0,1) N Q, let v € CE(X) be such 
that: 


(a) v(x) 2 « for x € B(x’); 

(b) v(x) =0 for x € X \ A. 
The countable set of all such functions is supposed to be a part of V, and hence 
they are taken into account in (5.14). Since each vs has compact support, for each 
compact À C X and any two configurations, ds(¥ N A, ¥ N A) > 0 only for finitely 


many s. Here 41 A := (p(¥) N A,n). 
For y ET, let a: y — A be a map, for which we denote 


la(z)) = M7 noo). (5.15) 


oca(z) 
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Then the pair 4 = (»y,a) is a marked configuration whose ground configuration is 
y and the mark map is a. By writing $ = (x,a) € ¥ we will mean that x € y and 
a € a(x). The configuration of marks a(x) = {a1,...Qa(x)|} yields the ages of the 
particles located at x € y, whereas |a(x)| is the total number of such particles. In 
some cases, we write a4 to indicate that a is defined on a given ^. Let Î denote 
the set of all marked configurations 4. Let also p : Î — P be the map such that 
D(y,a) = (y, |a), where |a|(x) = |a(x)| see (5.15). Then p o P maps 4 = (y,a) into 
its ground configuration y. For brevity, by writing p(^) we will mean (pop)(¥). Our 
aim now is to equip D with a complete metric. Define 


2 (reis (3, 7) 
KY, 4) = m l (5.16) 
- 1 ps Kan 4) 
Ks kn iQ 4) = | p NDS Wi s (o) — ` vs(x) ` viso). 
€p(* O€Q5 (a) z€p(4) ace 4! (a) 


Note that the latter can also be written as, cf. (5.13), 


(5.17) 


Js knl, a) = Vs(T)Wk n(o)- 
For a compact A C X, we write 4 1 A = (p(4) N A, a), where a is the restriction 
of a from p(¥) to p() NA. 


Proposition 5.1.6. For each € > 0, one may find a compact A, C X such that, for 
any two configurations, the following holds 


IKA, 4) — KAN Ae, 4 n Ae)| < e. (5.18) 


Proof. Fix € > 0 and then pick s, € N such that 2% > d Now let A, be covered by 
the supports of v, with s < są. For such s and all k,n € N, we have Ks knl, 5) = 
Ks kn (Y N Ae, 3 N At), see (5.16). This clearly yields (5.18). 


Since ø = 0 is in X, by (5.8) and (5.14) we have that 
d(B(), (5 )) € &(5./) (5.19) 
Proposition 5.1.7. The metric space (D, x) is complete. 
Proof. We begin by pointing out the following evident fact 


9— (s-ck-en),. ý (4 4^) 
R9, 3) :— > rr md d. 5.20 
T WR k,ncN 1 + Koenig s) i ut y) l ) 
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holding for all s € N and 4,4. Let now {4m = (Ym, am) fmen C Î be a &-Cauchy 
sequence. By (5.19) the sequence (D(^5) eu C Č converges to some ^. Take now 
x € p(y) and then pick a compact A C X such that AN p(y) = (xr). For this A, 
we then set 

Hore ` nyna ly), m € N. 


YEP(Ym)NA 


From the convergence of {P(4m)}men to ¥, it follows that nm(x) — n(x); hence, 
there exists m, € N such that n,,(x) = n(x) for all m > m,. Now we pick z' € X" 
and q € Q such that x € B«(z') and the closure of B,(x’) lies in A. Let now v, € V 
be such that v,(y) = < € (0,1) n Q, y € Ba(z’), and v,(y) = 0 for y € XX B(x’). 
For these m, and s, {4m}m>m,+1 İs also a &,-Cauchy sequence, see (5.20), for which 
we have 


sn (s osi) = 6 ^; MN Wado- P», $5 Wor nl@) 


YER(Am)NA GEA4, (y) YEP(Fm4 NA VEA4,, | (y) 


holding for all k,n € N, m > m, and l € N. 

Let us enumerate 2 = (r,o) € ĝm N A in such a way that apm X Qp+1,m for all 
p: This yields ĝm N A = ra Qim); so (Tnm, Qnm) } with n = huic) = n(x). 
Similarly, we have 


A ril NA= Trad O1,m41); DD (nmal, TAM 


with the same n. Then (01,5,...05,4] =: am € A, and also {a1 m+; Oma] =: 
am+1 € A, and the latter equality can be rewritten as follows 


TL 
Ks kn m, 270) E | px Vs, n (tpm) -y Uy n (Op m1) = SPk, nlQm, Am+1); 
p=1 


p=1 


see (5.12). By (5.20) and (5.12) we then get 


s 


plam, Gita) < sm. Aml): 


By Proposition 5.1.5 this yields the convergence of {am}m>m, to some a(x) € A, 
which holds for each x € ^. This defines the map a: p(y) — A, and hence 
the configuration 4 = (p(¥),a). Our aim now is to prove that k(4m,7) — 0 as 
m — +00. 

Fix £ > 0 and then pick a compact A; C X such that (5.18) holds with $ in 
the right-hand side. Let A? be its interior. Then pick compact mutually disjoint 
Az C A’, £ E p) NA? such that p(y) N A, = {x}. As A, is compact, p(y) N A? is 
finite. fiet (xj) jc; be an enumeration of it. For brevity, we will write A; in place 
of Az, j = 1,...,J. Similarly as above, by the convergence of [SC eR to ^, 
one finds My gud that p(y) NA, is a singleton and [D(45,) N A;| =: na, (v5) = n(zj), 


AT 


holding for all m > m, and j € J. Now we repeat the construction just made in 
each of A;. That is, we enumerate 


277 N A; e T Og) XE EN GE ceat 
and then set af, = C A Then we set 4am = (y. Axm), Where 
Yam = PI) NO A2 = (z1,...,2;) and asm(zj) = aj. In other words, the ground 
configuration of Ẹx m is the part of the limiting configurations p(Ẹ) contained in A?, 


whereas the marks are taken from the corresponding part of ĝm. By the triangle 
inequality we then have 


NET E A , at € 
K(Y, Vm) € KAN A2, Fam) + (4m N A2, Yam) + 3 (5.21) 


By (5.16), for each s, k, n € N, we have 


Kaien (FONE iam) = Sree) Y; vs) - wsale)]| 


L^ 
S 
WE 
“hd 
"7o 
co 
3 
Q 
3> 
= 
8 
S. 


Likewise, 


kinl OON om) — b» ICM) SEO 


I= Q€ y, 
J 
- Mus) Y, waala) 
j=1 acai, 
€ dim Ag YN Az) 


Both latter estimates yield 


KCA OAL, Jam) < J max p(aj, a(v;)) 
JS 


K(^m N A?, Aran) < Grr y). 


By the aforementioned convergence Jm — ^ and a2, > a(x;), one can find m, > m, 
such that the first two commands in (5.21) are smaller than ¢ for m > m,, which 
completes the proof. 
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5.2 Measures and functions on configuration spaces 
For v, and wz, as in (5.16) we set 
0. (1,0) = exp(—vs(2)wen(@)) — 1 = exp(—gs,an(#,@)) — 1, 


see (5.17). Then 0,4, (x, à) € C4(X) and 6,4, (x, 0) € (—1,0]. Let © be the subset 
of Ces( X) consisting of 


0(£) = e 9 — 1, = Ltn TUE ait (5.22) 
where the latter sum runs over a finite subset of N?. That is, each such g is a finite 


sum of gs kn defined in (5.17). Note that © is countable and closed under the map 
(0, 0") — (0 x 0"), where 


(0 x ONE) = 0(2) + 6'(2) + e(2)0 (2) = -14+ (14+ 0(x))(1 + 0'(z)). (5.23) 
Moreover, by (5.6), (5.8) and (5.22) it follows that 
g(x, o) € (v, 0), (5.24) 


holding for all a > 0 and x € X. Now for 0 € O, we set 


F*(4) = [[G +06) =e Y? (-s (8). Sef (5.25) 


LEY LEY 


Then F*(4) € (0,1] for all 4 € P, and hence F? € C,(D). The set of all such 
functions will be denoted by F?. For u € P(T), we then have 


it = [ Fou <1 


The Poisson measure 7, with intensity measure p satisfies 


T(E’) = exp (p(0)) = exp ( | &(£)p(d2)). (5.26) 


X 


For u1, u2 € PÔ) , their convolution is defined by the expression 
qu aX) = f Fé Uam (d) ua (53) (5.27) 
f 


that ought to hold for all F € B,(Î). For F° as in (5.25), it takes the form 


(H1 * uo) (F*) = i (F*)us (F^). (5.28) 
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Recall that a set F C C,(Î) is called convergence determining if ju, (F) > u(F), 
n — +00, implies jj, = u, holding for each (ji, eu C P(L). It is known, see Ethier 
and Kurtz (1986), [13], Theorem 4.5, page 113, that such F enjoys this property 
if it is closed under pointwise multiplication and is strongly separating. The latter 
means that, for each 4 € P and e > 0, there exists a finite family {F;} C such 
that 

inf max |Fj(7)— E3)| 20,  Ce=P\ BG). (5.29) 

4 €Ó. 
Here , 

BAY) = {7 : 854) < e, 

see (5.16). Note that taking € > 1 does not make sense as K(¥,7’) < 1 for all 
configurations. 


Proposition 5.2.1. The set F? is strongly separating and thus convergence deter- 
mining. 


Proof. By the very definition of O, cf. (5.23), F? is closed under pointwise multipli- 
cation. To prove (5.29), we note that, see (5.25) and (5.17), 


Pe gene (ans gam (Ren a (NP Y pice (neq) (5.30) 


holding for all 4,4’ € P. Now we fix 4 and e € (0,1) and then take m such that 
2-™ < $. For this m and any ^ € Ce, by (5.16) we readily conclude that 


A Af 
max Bx yo uk 
(s,k,n):stk+n<m 


which by (5.30) yields the proof. 


The important properties of the family F^ are summarized in the following state- 
ment. 


Proposition 5.2.2. The following is true: 
(i) BÊ) = e(F^); 
(ii) BD) is the bp-closure of the linear span of F*; 
(iii) F? is separating; 
(iv) F? is convergence determining. 


The proof of (i) and (ii) is standard, see Dawson (1993), Lemma 3.2.5 and 
Theorem 3.2.6, page 43, [11]. The proof of (iv) was done above, whereas (iii) is a 
direct consequence of (iv), cf. Ethier and Kurtz (1986), [13]. 


Proposition 5.2.3. The Kolmogorov operator introduced in (5.1) has the property 
L: F? 5 CL). 
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Proof. We consider each of the summands in (5.1) — denoted by L;,i = 1,2,3 sepa- 
rately. Thus, 


(LaF) = PG) | 6(0,0)x(de), (5.31) 
x 
i.e., it is just the multiplication operator by a 0-dependent constant. Next, 
(La F9)(5 - X im(2)0 £)F* (5X4) = F'(5)v5(4), (5.32) 
ic? 


where 


WA = 9 v) = 5 mae — 1). (5.33) 
Let us consider the following function 


o-(#) = g(£) — rm(#)(e%™ — 1), 7 > 0. 


Since each gs%n(%) < 1, see (5.17), it follows that g(2) < Je where Je is just the 
number of summands in the sum in (5.22). At the same time, m($) € m, « oo. 
Taking this into account, we set 


Then k 
br (£) > $ € X. 
Now by the simple inequality Be~7° < 4+, 7,8 > 0, we have, see (5.32) and (5.33), 
(LoF*)(4) < E | emeret (5.35) 
ic? 


which yields the boundedness in question. The continuity of LaF ? follows by the 
continuity of V5, which in turn follows by the fact that Y € C.,(X). Finally, 


(LFG) = (- Ya @) F6) = HOF). 


LEY 


ES vs; (x Jus, NG (o), 


which yields the continuity of Vj. At ds same time, by (5.9) we have 


By (5.22) we have 


(e, X à ser) wk, n,(@)| € õcg(z, a), (5.36) 


which yields 
lg (2, @)| € edo, (5.37) 


and also 


This completes the proof. 


We summarize the estimates obtained above in the following 
A ded. , 7C 

sup ((LF^)(5)| € x((,0)]) + m.e* ^ + —. 

4eb e 


Note that 
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(5.38) 


x(l&z, a) p= fi lé, a) |x (de) = f a-« (1-9) (dz) < L g(x, 0)x(dz), (5.39) 


see (5.24). 


5.3 The Kolmogorov equation 


5.3.1 Notions and useful estimates 


For 0 € O, see (5.22), we set 


6.(x, a) = 0(x, a + t) exp (Mta. a) — M(z,a- 2). 


Maa) =f m(x,6)d3 


(5.40) 


Let $ € X, then 2 > 0,(3) is continuous and compactly supported for all t > 0. 
Moreover, both maps t — 0,(ĉ) and a — 0,(x, o) are continuously differentiable and 


the following holds 


o 


o 
F —4;:(x, a) — m(x, a) Oz (x, a) 


ina) = 3a 


Note that 
(6,), (2) = 61, (8). 
Next, we define, cf. (5.22), 
g(%) = — log (1 + 6,(1)). 
By (5.40) it follows that 


l&(z,o)| < Gs, a +t)| =1—e 9) <1-e™%, 


where Jọ is the same as in (5.34). By (5.43) this yields 


dis). t>0, EX, 


(5.41) 


(5.42) 


(5.43) 


(5.44) 
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By (5.40) we also have 
[0 (x, o))| € 2m, + 0'(z, o + t)| € 2m, + Je, 
where m, is as in (5.34) and the estimate 
[9 (2, 0)| € |g' (x, 0)] < Jo 


was used, see (5.22) and (5.37). 
Now we define 


F?(4) = exp n (fa ae M@(da) ) do Fe (4, (5.45) 


with 0, as in (5.40). Clearly, F? € C,(Î) for all t > 0 and 0 € O, and 
0< F?) <1, Fel. (5.46) 


Furthermore, for all t, s > 0, the following holds, see (5.42), 


Fe, = exp gi i Bs 0)x(de)du) Fes. (5.47) 


Let us prove that F? € C,(Î) for all 0 € © and t > 0. As in the proof of Proposition 
5.2.3 we divide L into three parts. Similarly as in (5.31) we have 


(EFG) € x(06, 01) € x(aC,0)), 
see (5.39) and (5.46). Since g,(2) satisfies (5.44) for all t > 0, it follows that 
(LFA) € mee", 
holding for all t > 0, see (5.35). The estimate of |(L,F’)(4)| is obtained as follows. 


Denote oy 
q(x, o) = exp ( E f m(x, 8)d8). 
Then by (5.40) we have i 
$(g(r,a)) 2 G(x, a)®(g(z,a+t)), (b:-1-e*, b>0, (5.48) 


by which we get that g,(z,a) € g(a,a + t). Let us prove that 
gle, ae 9 > q(x, o)g(z, a + the 9, (5.49) 


By (5.48) this is equivalent to the fact that the function b > wA is decreasing, 
which is obviously the case. Now we take the a-derivative from both sides of (5.48) 


and obtain 


gilza) = (n o)g (o +t) exp (a(2,0) — gle, a 4 2) 


+ m(x,a + t)(exp (gir, o)) — 1) — m(x, o)(exp (gi(, o:)) — 1), 
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that can be estimated as follows 


ocg:(x,a)g'(x,a + t) exp (ac. a) — g(z,a + )) (5.50) 


gi, o?) 


+ 2m,(exp (gi(z,o)) — 1) € g(x, aje” (ac + 2m,), 
where we used (5.36), (5.44) and (5.49). Now we proceed as in obtaining (5.38), 
which eventually yields 
(LFG) < x(g,0)) + mee”! (2e 2m. = ly. —— (5.51) 


The key property of the latter estimate is that it is uniform in t. However, it does 
depend on 0. Along with the estimates derived above, we will use also the following. 
For 0 € ©, the corresponding g has the form as in (5.22). By (5.7) and (5.8) we 
have that 


V4 
exp ( — ~~ Jua (0) S Wise) < Wknl0), 
which means that, cf. (5.24), 
5V4 
egi. 0) gira) < g(a, 0). C:— exp (- 3 J (5.52) 


For ® as in (5.48), we have 
2 


which we use together with 5.52) to obtain the following 


n "Tom cg(x, 0 
a) > a)l) > exp(om.)e (TED) (833) 
eg(x, 0) eg(x, 0) 
Ia f Sa 
m EUM 
exp (—m.t)ceg(x, 0), 
and we have used the fact that g(z, 0) < Je, see (5.44). 


> exp (—m,t) 


IV 


where, Cg := m 


5.3.2 The operator 
We fix 4 € P and calculate the t-derivative of (5.41). This yields 


auno = (f $e vds) FG) + 3 Por NO (5.54) 


-: ‘ if 04(0,0)x(de)) EEG) + RA) 


icy 


do m4) [FPG \ 2) - FP @)] = (LEP). 


“LEY 


+ 
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This means that we have found a solution of the Kolmogorov equation for (5.1) in 
the following sense. It is a map t > F, € C (Î) , which is pointwise in 4 continuously 
t-differentiable and such that the equality in (2.4) holds. Our aim now is to solve 
(2.4) in a suitable Banach space. Recall that the paths t — 0, have the flow property 
(5.42), see also (5.47). Below, by saying of a property of 6,, s > 0, holding for all 0, 
we attribute this property to all 0, given in (5.40) with 0 taken from O. 


Proposition 5.3.1. For each 0 € O and s 7 0, it follows that Fb — F?s ast 0 
in the norm of Cy(T). 


Proof. For each 4, s > 0 and 0 € O, by (5.54) it follows that 


FG) - F6) = e (- [ f eso ote) FA.) - FPG (655) 


exp ( — / f bulz, 0)x(de)du) n / 7" (LS) da, 


which by (5.51) yields 


sup | Fe (4) — F%(4)| < tlo exp ( — [ l bu(2,0)x(dz)du). (5.56) 


4er 


This completes the proof. 
The next statement is a refinement of the one just proved. 


Proposition 5.3.2. For each 0 € O, s 7 0 it follows that LF” — LF” ast + 0 
in the norm of C,(T). 


Proof. First of all we note that the equality in the first line of (5.55) allows one to 
obtain the property in question by showing that LF?,, — LF? in the same sense. 
As in the proof of Proposition 5.2.3, we split L into three parts and consider each 
of them separately. 

Fix 0 € O and then denote 


itu) = [sor MO = f itat (6.87) 
Then 
file + sent — eo) 
lie" E, (4) — F9) 
=: I,(t)+ I(t), 


(aE I — G8 F2)(9)] 


IA 


+ 


see (5.46). Then A(t) — 0 as t > 0 since 7(t)e" is a continuous function of t. At 
the same time, J9(t) can be estimated as in (5.56). This yields the proof for L3. 
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Now we proceed to La, for which it follows that 


(L2Fp")(7) — (25 F^)(5) = JG) + Kâ), (5.58) 


OUR A E 


By (5.45) and (5.57) we have 


HG UG), sen qoe c qam 6n (5.59) 
v oy) e aps YQ) ee). 
By (5.39) and (5.40), (5.57) it follows that 
jeher nes Te motos 0); 
which then yields 
Tilt, P¥) < tx(9(-,0)) exp ( — 9 gere (ê)) (5.60) 


To estimate Y» we write 
h, (2) = min {g+ (2); g.(1)). 


Then by (5.50) and (5.53) we get 


Talal € exp (— So hela) lass) — 0G) (5.61) 
< exp (Joy (2c + 2m.) exp ( — 2 hel (8)) )»E[ u(x, a)du 
< exp (Jo) (Tc + 2m.) exp ( upto ` g(x, 0)) »3 tg(x,0) 


LEY BEY 


= tC 0 0 Yo qy (4), 


56 


At the same time, 


^ 


—-m(2)0... (1) J J ^ 
0< <m,e°? A =: me ?WVi(t, ^). 
< Dee is Dek &) » d i( 4) 


Thereafter, we have 
|J(4)]. S tre” xo 0) Mune 99 
+ tCom,e” Y(t, 4) Vo(4)e POM 
= h(@,7) +b, 4), 
where C», Co(t) and Wo are as in (5.61). Then 
II, (t, 7) € tme !x(g(.,0)) — 0, t + +00. 


Let tọ > 0 be the (unique) solution of 


e "ta, = 94/t. 
Then for tg, we have 
Ib(t4) € tCsm.e^[Wo(4)]^e (099.0) (5.62) 
« WtCom,e7~? exp ( — (e(t) — 2/09) 
€ viComue ^, 


which yields the convergence 


sup |J;(4)] — 0, t Q0. 
ser 


Now we turn to K,(4). First, by (5.41) we have 
E — Or40(2)| = | ime Em ma, a)&,(z, o))du|. (5.63) 
Next, by (5.43) and (5.40) it follows that 


eH) |g (a, o)] 


(x ,a) 


p 


IA 


e^ (Ge + 9m,)gu(x, a) 


< e” (oc ein 2m, )g(z, 0), 


and also 
|m(x, a)8,(z, a)| € m.g«(v, a) € m.g(z, 0). 
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The latter two estimates yield 


LHR(5.63) < e” (€+ 3m,)g(z, 0). 


By (5.52) this yields 


KU) 


tm,(8 + 3m,)e" wo (4)e- 9*9? 


IA 


e3Je-1 


— 0,t — +00. 


IA 


tm, (ao + 3m,) 


By (5.58) this completes the proof for L3. Next, we write 


Then 


(Li Fy*)(4) — (Li F"6s)(4) = QA) + Rel), (5.64) 


+ F%(4) S 0nd) [oA ES 1 


=: QP) +P A). 


By (5.50) and then by (5.59), (5.60), (5.61) we get 


RPA) 


Then 


< 


IA 


e^ (5c + 2m,)(Ys(t, 4) + Talt, 4) 2 gus (ê) 


LEY 


and also, cf. (5.62), 


=P (3) 


IA 


VtCse"* ? (ac + 2m,) exp ( (C9 (t) = 2/4) 


< VtCsm,e" ? (ac 4 2m,), 
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for t € tg. The latter two estimates yield 


supiQ(4 20, tO. (5.65) 
ser 


Next, by (5.44), (5.50) and (5.52) we have 


QP AL € AFA) Y dirlo) — 9s (6) (5.66) 
LEY 
t+s 
< eMe (oc-- 2m,)°F" (4 Op» gis 22] g,($)du 
LEY d 
< tJa” (ac + 2m,)?Wo(4) exp (—@Vo(4)) 
z 2 
« m anti d — 0, 0: 
[e 
which together with (5.65) yields 
sup |Q«(5)| 0, t 0. (5.67) 


4ef 


Now we turn to estimating R;. By (5.40) we have 


1.,(2, a) -6(,0)| < |O(t,a+t+8)||G+s(v, a) — G(x, @)| (5.68) 
+ j0'(z,o c t4 s) — 0'(r, o — s)| 
*o[m(zx,a-ct-s)-—m(xz,o-s)|-]|0:.5(x,0)]| 
*om(z,os)-m(z,a)[0.(v, a) 


=: 61(t, £) + do(t, 2) + d3(t, 2) + 04(t, £). 
By (5.50) we have 
a+t+s 
(2) € em pg ates f m(e,8)d3 (569) 


< tm,e" (ac + 2m,)g(x, 0). 
To estimate à», by (5.43) we first get 


e(&) = -g (@)e™, 
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by which we then obtain 


6x(t,@) < ]|g'(zx,a - t - s) — g'(z, o + s)|e Iet) 
«gs, + s)||e 9 (5o***9 — e7s (mes) 
< |g(z,a+t+s)-—g'(z,a+ s)| 
+ [g(r,o- s)||g(z,a +t +s) — g(x,a + s)| 


= 05.1(t, i) + 0 2(t, m). 


Now we recall that g(z,o) is as in (5.22) with Wk, defined in (5.8). Thus, we 
can write 


t+s 
bati) You, (o) [jw "hy, nio eu. (5.70) 
j S 
For each k and n, we have 
(a)e- 7*9) + [opun (a)]^e t| 


ou, (o)| + |a, (a?| € C, 


^ 
A 
g 


holding for some C > 0 that is independent of k,n and a. The latter conclusion 
follows by (5.9) and the fact that |u, (o)| < 26(no?) with 


EMIT. 
o(B) = a+ 8y B 2 0. 
Then by (5.70) we get 
ôz1(t,ĉ) < tCg(e, 0). (5.68) 


At the same time, by (5.36) and (5.37) it follows that 


t+s 
alia) e doors] f lg'(w, a+ u)|du < t(Ge)?Jeg(x,0), 


which together with (5.3.2) yields 
ó»(t, 2) < t[C + (ac)? Jolg(z, 0). (5.68) 
Finally, (5.2) and (5.43), (3.13) we have 
b3(t,@) < x(t) ge(£) € x(t) g(a, 0). 


The same estimate holds true also for d4(t, #). Then by (5.68) and (5.69), (5.3.2) we 
have that 
(6: (7,0) — O'(x,a)| < w(t)g(w,0), w(t) 30, t0. 
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holding for some continuous function w and all @ € X. Now we use this in (5.64) 
and obtain, cf.(5.66) 


(RA) S w(t)” F (A)T) € we” Wo (4)e MM < 


l (5.68) 


which together with (5.67) yields 


sup |(LiF”)(7) — (LiF*)(7)| ^ 0,t — 0. 
4cr 


This completes the whole proof. 


5.3.3 The domain 


We recall that © consists of the functions as in (5.22) and the countable collection 
F? C C,(D) consists of the functions introduced in (5.25). It has a number of useful 
properties established in Propositions 5.2.1 and 5.2.2. Let Co be the linear span with 
rational coefficients of the set F^ : s € Q,,0 € O, i.e., each F € Cp is a finite linear 
combination of F^, with positive rational s and 0, given in (5.40) with all possible 
choices of 0 € O. 


Remark 5.3.3. The set Co is countable. It enjoys all the properties mentioned in 
Proposition 5.2.2. 


Now we set 
C=; (5.69) 


i.e., C is the closure of Co in the norm of Cy(D), which we denote || - |. With this 
norm it is then a separable Banach space. 
For A > 0 and 0 € © and s > 0, we define, cf. (5.45), 


Elta f eX FP-(4)dt (5.70) 
0 


a exp |- At + n ( f 0ste, aje E ta))do] Fo (4)dt. 


Since 0,(x, o) < 0 and F?* satisfies (5.46), the above integral converges for each 4. 
By the dominated convergence theorem and the boundedness Fy" (4) < 1 it follows 
that Pis, € C,(Î). Moreover, it can also be understood as the Bochner integral 
in the latter Banach space. Therefore, Fà can be approximated in ||- || by the 
Riemann integral sums centered at rational t, which means that 


Fy 9, € C, for all s > 0 and A > 0. (5.71) 


At the same time, we also have 


1 i 
0 < Frye. lã) < Y VET. (5.72) 
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Let o € (0, +00). The continuous differentiability of the map 


A 


a — Fyo(V¥\ ĉU (a,a)) € 


for each ^ and ĉ € ^ follows from the dominated convergence theorem. 
Lemma 5.3.4. For each À > 0,s > 0 and 0 € O, the following holds 

LFyo, = AFy, — F^ (5.73) 
Proof. By (5.45) and (5.70) we have 


-Foo -Foo 
LFye, = L | edt f e LF” dt (5.74) 
0 0 
+00 ð 
= f e Fbdt = —F% + AF, 
ü Ot 


where we have taken into account the upper bound in (5.72). The commutation 
L f = f L can be justified by means of the Lebesgue dominated convergence theo- 
rem. 


Lemma 5.3.5. For each 0 € © and s > 0, it follows that |XFy o, — F% 
[ALF o, — LF%|| > 0 as A — +00. 


Proof. In view of (5.72), {AF}, : A > 0} is bounded. By (5.70) we have 
+00 et 
AP.) = f exp (—t4 f Cf tasela0)x(dz)) da) rnt 
0 0 X 


where £ := A— 1. 
Then by (5.56) it follows that 


— 0 and 


Too l 
bb - Fas [^ ek) - FO 
0 


which yields that ||AF\,9, — F% 
5.3.2 we have, cf. (5.3.2), that 


— 0 as A > +00. In the same way, by Proposition 


1 
JALFyo, — LF*|| < a(5) zi; X. eq, 


holding for an appropriate continuous à such that &(e) > 0 as € > 0. 


Let Do(L) denote the linear span of the set {F}, : À > 0,s > 0,0 € O}. 
By Lemma 5.3.5 and (5.69) it follows that 


Ese Did), (5.75) 


i.e., Co is contained in the closure of Do(L) in the norm of C;,(L). 
Hence, Do(L) is a dense subset of the Banach space C, see (5.69). Define 


FIL = FI (ZF I], F € DL) (5.76) 


where as above ||- || is the norm of C,(I) . 
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Definition 5.3.6. By the domain of the Kolmogorov operator L, denoted by D(L), 
we mean the closure of Do(L) in the graph-norm introduced in (5.76). 


Lemma 5.3.7. The operator (L,D(L)) is closed and densely defined in C. Its re- 
solvent set contains (0,+00) and Co C D(L). 


Proof. In view of (5.71) and (5.75), the ||- ||-closure of D(L) is C. The closedness 
of (L, D(L)) is immediate, and the inclusion Co C D(L) follows by the second part 
of Lemma 5.3.5. By (5.73) it follows that the resolvent of L, denoted R)(L), has 
the property 

RA(L)F* = Fio, A>0, s>0, 0€86, 


by which and (5.72) we also have that the operator norm of R)(L) satisfies 


1 
IROI S 5 


as F?: form a dense subset of C. This completes the whole proof. 


5.3.4 Solving the Kolmogorov equation. 


The result obtained in Lemma 5.3.7 allows one to solve the Kolmogorov equation 
(2.4) in the following sense. 


Theorem 5.3.8. Let (L,D(L)) and C be as in Lemma 5.3.7. Then, for each F € 
D(L), there exists a unique continuously differentiable map 


[0, +00) 251 — F; € D(L) CC, 


which solves (2.4) with Fy = F. In particular, for F = F**, with F® as in (5.40), 
s 2 0 and 0 € O, the solution has the following explicit form, cf. (5.45) 


FG) = exp ( | l ( a 6, (1. o) x(dz))da) F** (4). (5.77) 


Proof. By the celebrated Hille-Yosida theorem, see, e.g, Pazy (1983), page 8, [32], 
(L, D(L)) is the generator of a Cy—semigroup of bounded linear operators S(t) : C > 
C such that ||S(t)|| = 1 and the solution in question is F, = S(t)F, the uniqueness 
of which is also a standard fact, see Pazy (1983), Theorem 1.3, page 102, [32]. The 
validity of (5.77) follows by the calculations as in (5.54). 


5.4 The result 


5.4.1 The martingale problem 


We begin by recalling, see Proposition 5.2.2, that the class of functions Fo, see 
(5.25), is separating, i.e., if j4(F) = we(F) for all F € Fo, then jj = u2, that holds 


63 


for each pair u1, ji» € P(D). Next, for t > 0, € P(P) and 0 € ©, we determine 
ut € P(T) like in definition 4.1.2 by the relation 


p (F°) = u(F*). (5.78) 


Recall that, for a positive Radon measure 9 on X, the Poisson measure 7 with 
intensity measure ọ is defined in (5.26). For t > 0, we then introduce the Poisson 
measure Ti = To, by defining its intensity measure 


0:(dz) = 1o4)(a) exp (C M(£))x(dx)da, (5.79) 


where x is the same as in (5.1), da is the Lebesgue measure on R+, M is as in (5.40), 
and 1,9 4)(-) in (2.1). Note that mo({Ø}) = 1 since z9(F?) = 1 for all 6. Then we set, 
see (5.27), 


Jg = Tik pl, t>0. (5.80) 


Let à4 be the Dirac measure centered at a given ^ € Î. Then 
pi: mex 05 (5.81) 


is a transition function, cf. Ethier and Kurtz (1986), page 156, [13]. Indeed, pi € 
P(L), pj = ð, and the measurability of the map (5,4) > p/(B) € R, B € B(T) 
, follows by the measurability of (t, 5) + 64(B) € R and the continuity of t > 
7,(B) € R. In view of the separating property of F°, see Proposition 5.2.2, item 
(iii), the flow property of {pi asi can be obtained by showing that 


p(y fE w ts>0, OEO. (5.82) 
T 


By (5.27) we have 
rite) = es (f. f tsi) 
exp TL balz, 0)darx (dee) ) F^ (4^), 


see also (5.78) and (5.40). By the latter formula and (5.79), (5.81) we then get 


RHS(5.82) = exp ( f f PR CAES ? f ‘Prale, ydo (dr) 8s CP) 


— exp (f. a ba(z, 0)dax (der) ) F^ (4) = LHS(5.82). 
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As is known, cf. Ethier and Kurtz (1986), Theorem 1.2, page 157 [13], the tran- 
sition function (5.81) determines a Markov process A with values in I, the finite- 
dimensional distributions of which are given by the following formula 


P(X(5) € By,..., (5, -] Js 7 pict (B uL (ds) 
Bı Bn-1 


x X2 au (dT2)p), (d51)u(d3), (5.83) 


holding for all n € N,0 < s; € s3 € ... € s, and B; € B(Î). Here u € P(D) is the 
initial distribution of X. Our aim is to show that such a process is unique up to 
modifications. 


5.4.2 The statement 


The process determined by (5.83) describes the stochastic evolution of the population 
which we consider. To verify whether it is the only one, we have to specify which 
processes of this kind can be associated to the model defined by the Kolmogorov 
operator (5.1). As is standard, the corresponding specification is made by their 
martingale property, see Ethier and Kurtz (1986), Chapt. 4, [13]. 


Definition 5.4.1. Let X be a measurable process on some probability space (Q, S, P) 
with values in T. Let also {§:}is0 be a filtration such that X(t) and 


/  G(X (u))du 


are $,— measurable for all t and G € BÔ). We say that X is a solution of the 
martingale problem for (L,D(L)) if for each F € D(L), 


M(t) = FXO) - | LPW) 


is a &1— martingale. If there exists a solution of the martingale problem for (L,D(L)) 
and uniqueness holds, we say that the problem is well-posed. In the same way, we 
define the martingale problem for (L, D(L), u) if the initial distribution u € P(T) is 
specified. 


The process related to the transition function (5.81) solves the martingale prob- 
lem for (L, D(L)). Its uniqueness will be shown by proving that all other solutions 
have the same finite-dimensional marginals, i.e., they coincide with those defined in 
(5.83). We are going also to show that the solution is temporarily ergodic. 


Definition 5.4.2. Let the martingale problem for (L, D(L), u) be well-posed. Then 

LE T(l ) is said to be a stationary distribution if for each n and 0 < sı < s3 < 

. € Sn, the n—dimensional marginals introduced in (5.83) corresponding to t + 
.,t+ Sn are independent of t > 0. 
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Definition 5.4.3. By a solution of (2.3) we understand a map R4 Ə t > m € P(L) 
possessing the following properties: 


(a) for each F € C,(ÛÎ, the map Ry > t > ju(F) € R is measurable; 
(b) the equality in (2.8) holds for all F € D(L). 


If P is as in (5.83), then u is stationary if and only if 
u= i pip(d4), 
f 


holding for all t > 0. Now we can formulate our result. 


Theorem 5.4.4. The martingale problem for (L,D(L)) is well-posed in the sense of 
Definition 5.4.1. Its solution is defined by finite-dimensional marginals, see (5.83), 
with the transition function defined in (5.81). If the departure function satisfies 
m($) > mo > 0, holding for all $ and some mo, then there exists a unique stationary 
distribution u = m$, which is the Poisson measure with intensity measure 


o(dz) = exp (-M(2))x(dx)da. (5.84) 


Moreover, in this case the solution A of the martingale problem for (L, D(L), u) is 
temporarily ergodic in the following sense. Let u € P(T) be the law of X(t), t > 0. 
Then ut = To as t — +00. 


Proof. The proof of Theorem 5.4.4 is divided into the following steps: 
(a) Proving uniqueness. 
(b) Showing the stationarity and ergodicity if m(z) > m, > 0. 


The realization of (a) is based on the Fokker-Planck equation (2.3) for this for L. It 
turns out that its solutions can be obtained explicitly. We are going to do this now. 


Lemma 5.4.5. For each uo € P(D), the map t — u defined in (5.80) is a unique 
solution of (2.3). 


Proof. For each 0 € ©, the map t > uu (F?) = ug(F**) is continuous by the domi- 
nated convergence theorem and hence measurable. By (5.28) we then have 


pF?) = mF?) pe F%) = mF?) uo (F*) 


exp ([ | f 0 07 az) da) nol”); 


Thus, the map t > pu ( F°) is continuous and hence measurable. Then the measura- 
bility of t > iu (F) for all F € © (T) follows by claim (ii) of Proposition 5.2.2. Now 
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we turn to proving the equality in (2.3) for F = F4,5,0 € O, see Definition5.3.6. By 
(5.80) and (5.45) for s,t > 0 and 0 € O, we have 


ns) = exo (f ( f bale Oxido) da) aE trt 

ee Ti ( [ 6r. 0)xlda))da 

f Cf tte) da’) wo(#") 

exp (f (f bate.0)xde)) da) (r^) = n (rt... 


Then by Fubini's theorem and the latter fact we get 


+ 


nO) = Hol Fro) = fe [eco - nol F2)] a 
0 
en Fe ,— F? ds 
li of t+ ) 
+00 t ð 
—As fidt F? 
li e / zrol 2, )dsdu 
+00 t ð 
B ud d — hu F? dsdu 
f pm 
-Foo 


' as 9 pe i 
= nu e 5 Fas) - f Hul LF, o)du, 


where we have used also (5.74). Now we prove uniqueness by applying arguments 
similar to those used in Costantini and Kurtz (2015), Lemma 2.11. Assume that a 
map t — Ju satisfies condition (a) Definition 5.4.3 and F, G € C(I) are such that 


t 


WE aR) = J E 


holding for all t > 0. Then the map t — ju(F) is almost everywhere differentiable 
and 
dp (F) = ui (G)dt. 


Then integrating by parts we get 


2 t 
E | e™y,(F)ds = e-™uy(F) = uo(F) — " e p. (G)ds, 
0 0 
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which yields 


lh) = eu) f e DasF) — sal) ]d, 


holding for all £, À > 0. Passing here to the limit t — +00, for F = Fs and 
G = LF) o, see (2.3), we arrive at 


coo Teo 
uol F9) = f e "u,(AFy — LFy9)ds = f e ^*y,(F^)ds, (5.85) 
0 0 


see (5.73). Assume now that (2.3) has two solutions, u; and ft, satisfying the same 
initial condition ju|, y = ft|,.g = Ho. By (5.85) the Laplace transforms of both maps 
t > p,(F°) and t > (F°) coincide, which yields (F°) = (F°) holding for each 
t and all F?, 0 € O. Then the uniqueness in question follows by Proposition 5.2.2. 
This completes the whole proof. 


The existence of a solution of the martingale problem for (L, D(L)) was shown 
by the very construction of the finite-dimensional marginals of X in (5.83). To prove 
uniqueness we use the following fact, see Ethier and Kurtz (1986), Proposition 4.2, 
page 184. Given u € P(D), let ¥ and X' be solutions of the martingale problem 
for (L, D(L), u) whose onedimensional marginals, 4, and ji^, coincide for all t > 0. 
Then all their finite-dimensional marginals coincide and hence the problem is well- 
posed. Clearly, both u; and u, solve the Fokker-Planck equation with the initial 
condition u. Then they coincide by Lemma 5.4.5. This yields well-posedness. 

Now we show the stated ergodicity. If m satisfies m(ĉ) > mo > 0, then M(x, a) > 
moa, see (5.40), which for 0 as in (5.40) yields 


VesolOe(#)| < e 9. 


By the continuity of the map t > u(F**) we then get 


P) = (f [esos ox tz) ta) uF) (5.86) 


— exp QU [ f 6s 0)x(de)]da), t + +00. 


By claim (iv) of Proposition 5.2.2 this yields u; = 7, see (5.84), holding for each 
initial u € P(L). Clearly, uj = 7, if i = To, which means that 7, is a stationary 
state. If there exists another stationary state, say u’, then (5.86) fails to hold for 
u = u', which contradicts the convergence just established. This completes the proof 
of Theorem 5.4.4. 
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Chapter 6 
Appendix 


Here we prove that 7 introduced in (5.3) satisfies (5.4). During the whole proof, we 
deal with the function l(a) = a+ i for a > 0. In case ag = 0, and 04,05 < 1 we 
have w(a,) = 04, w(a2) = as, r(o1,03) = |o — a,|. The triangle inequality turns 
into true statement: 

ay € |as — o1| + a3. 


For a3 > 1 we can have: 


1 
Og — Q1 < Qı + — 
Q3 


or 


1 
Q3 — Qı > Qı —. 
Q3 


For the first one we get: 


1 
wai) =a, < az — a + B = 7r(ay, a3) + w(as). 
3 


Then 
204 < I(a3). 
The equality holds only for o; = a3 = 1. In the second case we have true inequality 
o1 € a + m 
For o4 > 1, a3 € 1, we have two possibilities: 


1 
a1 — Q3 < &3 + — 

Qı 
or 


Q1 — Q3 > Q3 + —. 
Q1 


In the first case we have 
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which yields (5.4). In the second case we have: 


1 1 
wla) = — <a34 F Q3. 
Oy Oy 


It remained to consider the situation when o3 > o3 > 1 and 1 < o4 < az. 
When qa; > a3 > 1 for ay — a3 < acti we have: 


which proves (5.4) for this case. 
When a3 > a3 > 1 for a; — a3 > me we have: 


which yields (5.4). 
When 1 < o4 < a3 for o4 — o3 € zn + a we have to prove : 
1 
zc ugue 
Q1 Q3 
which is equivalent to I(a,) < l(o3). The latter follows by ag > a, > laslis 
increasing. For a3 — a, > m + ds the proof of (5.4) is immediate. 
Now we consider the case 0 < o4 < o». If ag € 1 and o3 = 0, then 


Ay — A, € w(ay) + wla) = o1 + oa. 


For o5 € (0, 1], it follows that 7(a3, a;) = |a3—a;|. Then (5.4) turns into the triangle 
inequality for | - |. 
The same is true also for o3 > 1 such that ag — F < 2a,. For 


204, < o3 — ES < 202, 
Q3 
the right-hand side of (5.4) is ag — ag + o4 + oe Then az — o4 € RHS(5.4) turns 
into 2(a2 — o4) € l(aa), which holds since 2(a@2 — o1) € 2 < l(aa) for a3 > 1. For 
205 < 03 — 25 the right-hand side of (5.4) is o + ag + 2/a3, which is bigger than 
Q2 — Qı. 
Consider now 0 < o4 € 1 < o» and ag — al < 2a,. The latter means that 


f (ao, 01) = Qə — Qı. 


For a3 = 0, the right-hand side of (5.4) is o4 + 2 and the latter turns into o» — 
= < 2a,, which holds in this case. The same is true for 203 < ag — m For 
Q2 — a < 2a3 < 201, the right-hand side of (5.4) is ag — a3 + o4 — a3, which is 
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bigger than o» — o4 as a3 € a4. 
For o4 € a3 < 1, the right-hand side of (5.4) is 


Q2 — Q3 + Q3 — Q1 = LHS(5.4), 


Next, consider 1 < o3 < a», where 


1 1 
Og — — Oo — — < 2a. 
Q3 Q2 


For ay — = < I(a3), the right-hand side of (5.4) is 
Q2 — Q3 + Q3 — Q1 = RHS(5.4). 


The case of a2 — T > I(a3) > 2 is impossible since ay — em < 2a, € 2. For a3 > ag 
such that a3 — + < 2a, and o3 — + X l(o3), we have F(a3,a,) = a3 — a, and 
7 (a3, @2) = a3 — a». Then the right-hand side of (5.4) is 


a3 — Q1 + a3 — Q2 = 203 — (a2 + 01), 
which is bigger than o» — o4. The case of a3 — m > l(a3) is impossible for a3 — 
a < 2a,. For a3 > ag such that a3 — F > 2a, and a3 — em < l(a3), we have 
F(a3, 04) = 04 + ad and 7(a3, 093) = a3 — o». Then the right-hand side of (2.2) is 


I(a3) — a2 + O1, 


which yields (5.4) in the form 2(a2—a1) € l(a3). By o3 — al < 205, a2 > 1, we have 
that 2(a2—a) < 24/1 + aj, whereas a3— - > 2a; yields l(a) > l(o.) = 2/1 + aj, 
which proves (5.4) in this case. Here a, is the positive solution of a — i = 2a. 
For a3 > a» such that dcc > 2a, and ag ce > [(a3), we have (a3, a1) = ages 
and T(a3, 05) = T + a Then (5.4) turns to 


1 
az — a, € ay + 2/az+—, 
Q2 
which holds as ag — al < 2a}. 

Consider now 0 < a; € 1 < o» and o» — al > 2o,. The latter means that 
f(01,02) = ats. For a3 € [0, a1], we have that F(a, a3) = a1—agz and F(a2, a3) = 
a3 + on Hence, (5.4) turns into equality. 

For a3 € (a4, 1] such that a; — al > 203, we have the right-hand side of (5.4) in the 
following form 

1 

Q3 — 01 + Q3 + os 

a2 
which is bigger that 7(04,02) since o, < a3. 
For o5 € (a1, 1| such that ag — = < 2a3, we have that the right-hand side of (5.4) 
is 

a3 — Q1 + Q2 — a3 = Q2 — Q1 > f(01, 02). 
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Consider now a3 > 1 such that a3 — al < 2a,, which means that a3 < ag 
and f(o4,03) = a3 — o4. For F(a2,a3) = a» — oa, the right-hand side of (5.4) is 
Q2 — a3 + a3 — 04 > F(A2, a3). For 7(a2,a3) = us + a the right-hand side of (5.4) 
is 


hence, (5.4) turns into 204 < l(a3), which holds since o4 € 1 and l(o3) > 2. For 
a3 > 1 such that a3 — ds > 204, we have that (04,03) = o1 + T Then (5.4) turns 
into 

1 NER 
— X — + f (Qe, Q3); (6.1) 


a2 Q3 
which clearly holds for a3 € o, and also for a3 > o», where for (05,03) = a3 — a» 
it turns into /(a3) < l(a3) — which is true as l(a) is increasing for a > 1. For 
F(a, 03) = E + a) the validity (6.1) is immediate. 

Let us consider now the case of 1 < a, < @ and l(o1) € ag — ae where 
f(a5,01) = m + oat For a3 € 1 such that 204 < a, — = we have 7(a3,Q;) = 
a3 + ay = 1,2. Then (5.4) obviously holds. 

For a3 < 1 satisfying a, — a < 203 < 05 — € we have 7(a3,04) = o1 — a3 and 
F(a3,Q2) = a3 + a Hence, (5.4) turns into equality in this case. The remaining 
case Qs — a < 2053 < 2 is impossible since o» — = > l(a) > 2. 
For a3 > 1 such that 2a3 < ay — T we have 7(a3, aj) = al + + i = 1,2. Then (5.4) 
obviously holds. For a3 € (1, a1] satisfying o1 — m < 203 < a2 — A (5.4) turns 
into 
1 1 1 1 
— F— <o,-03+—+4+-, 
Q1 a2 Q2 Q3 
which holds for ag < o4 as the function a — i is increasing. For a3 € (a1, a2] such 
that Q3— a. < l(a) < Llas) < Q2— =, the right-hand side of (5.4) is Cot sara T 
and hence the latter turns into I(a,) < l(a3). For a3 € (a1, Q2| satisfying a3 — oa < 
l(o4) and ag — al < l(a3), we have (5.4) in the form al 4 m X ay — Q3 + Q3 — 0, 
which holds as (œi) < o» — aed For a3 € (01,03] satisfying l(o1) < o3 — oe and 
a2 — al < I(o3), the right-hand side of (5.4) is 
1 1 
— + — +2 — 03 
Qı Q3 
which is bigger than 7(a1, o3) as a—+ is increasing. For a3 > o» such that a3 — as < 
I(@2), the right-hand side of (5.4) is 
1 1 


— + — +03 — a3. 
O1 Q3 


Hence, (5.4) holds as l(o3) < l(az). For a3 — em > [(a3), (5.4) turns into 
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Finally, let us consider the case of 1 < o4, < o» and l(o1) > ag — where 
F(Q2, 1) = ag — o4. For a3 € 1 such that 2a3 < 01 — m it follows that T(o;, a3) = 
as + 1, i — 1,2. Then (5.4) turns into 


1 
Qo. — — < l(a) + 203, 
Q2 
which evidently holds in this case. For o4 — a < 203 < ag — m we have that 
f(01,03) = a, — a3 and 7(a2,a3) = a3 + i Then (5.4) turns into 


1 
QU a < 201, 
Q2 


which is the case for 204 > l(04) > o» — A 


For o4 € 1 satisfying 204 > o» — = it follows that 7(a;,a3) = a; — 03,0 = 1,2. 
Then (5.4) turns into 

203 < 204. 
which is obviously the case. 
For a3 € (1, 01] such that /(o3) € ay — m. we have that ?(a;, a3) = +,+ ad = 12 
Then (5.4) amounts to az — ES < I(a,) + 2/a3, which obviously holds. 
For a, — m < I(a3) € ag — = we have that 7(a1,a3) = o1 — a3 and 7(a2, a3) = 
= + ae Then (5.4) amounts to 

1 


1 
ai — Qas Oi + — 2 a- —, 
Q3 Q2 


which is the case for 
1 1 1 
a, — Q3 + a1 + — > œ + — > l(a) > Q — —. 
Q3 Q3 Q2 
For az — al < l(o3), we have f(o;, 03) = a; — 03, i = 1,2. Then (5.4) is 


a2 — 05/55 Q2 — Q3 + Q1 = Q3, 


which obviously holds as a3 € o1. 

Now we consider a3 € (a1, a5]. For as — + € l(o1), we have (a1, a3) = a3 — o4 
and T(a9, 03) = a» — a3, which yields equality in (5.4). 
Recall that ay — m < l(a); hence, a3 — 7 > [(o4) is impossible for a3 < o5. 

It remains to consider a3 > a2. For az — oA < I(a,), we have that F(a;,a3) = 
a3 — a;, i = 1,2. Then (5.4) takes the form 


a2 — Q1 = 203 — Q1 — Q2, 


which obviously holds in this case. For l(a,) < a3 — al < I(a3), (5.4) amounts to 
203 < I(a3) + l(a), which holds as l(o1) > o» — + (assumed) and l(o3) > a2 + a 


a2 
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for a3 > a». For a3 — al > (az), we have that ?(03,0;) = 


I-1,i—1,2. Then 
(5.4) turns into 


c3 


1 1 
a9 — a, € 2/03 + — + —, 
a2 [941 


which holds since /(o1) > a2 — +. 


This completes the whole proof of Proposition 5.1.1. 
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